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ON THE REGULARITY FOR THE NAVIER-SLIP THIN-FILM EQUATION IN 

THE PERFECT WETTING REGIME 


MANUEL V. GNANN 

Abstract. We investigate perturbations of traveling-wave solutions to a thin-film equation with 
quadratic mobility and a zero contact angle at the triple junction, where the three phases liquid, gas, 
and solid meet. This equation can be obtained in lubrication approximation from the Navier-Stokes 
system of a liquid droplet with a Navier-slip condition at the substrate. Existence and uniqueness 
have been established by the author together with Giacomelli, Kniipfer, and Otto in previous work. 
As solutions are generically non-smooth, the approach relied on suitably subtracting the leading- 
order singular expansion at the free boundary. 

In the present note, we substantially improve this result by showing the regularizing effect of the 
degenerate-parabolic equation to arbitrary orders of the singular expansion. In comparison to related 
previous work, our method does not require additional compatibility assumptions on the initial data. 
The result turns out to be natural in view of the properties of the source-type self-similar profile. 
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1. Introduction 

1.1. The thin-film equation as a classical free-boundary problem. We are interested in the 
thin-hlm equation with quadratic mobility 

dth + dz =0 for t > 0 and 2 ; > 2’o(t). (1.1a) 

This is a fourth-order degenerate-parabolic equation modeling the evolution of the height h of a 
two-dimensional thin viscous him on a one-dimensional hat substrate as a function of time t > 0 
and base point z > Zo(t) [6,8,27]. For simplicity we assume that the droplet extends inhnitely 
to positive z and has a free boundary 2’o(t) denoting the contact line, that is, the triple junction 
between the three phases liquid, gas, and solid (cf. Fig. 1). Necessarily 



h = 0 for t > 0 and z = ZQ(t), (1.1b) 

which simply sets the location of the free boundary to be Zoit). Secondly, we assume 

dzh = 0 for t > 0 and 2 ; = Zo{t), (1-lc) 

leading to a zero contact angle at the triple junction, known as complete (or perfect) wetting. The 
notion can be explained if one considers quasi-static droplet motion in which the contact angle is 
determined by a balance between the surface tensions of the three interfaces at the contact line 
{Young’s law). The condition dzh = 0 at z = Zo{t) implies that an equilibrium is generically 
not achieved and therefore the him will ultimately cover the whole surface. Finally, a condition 
determining the evolution of Zo{t) is imposed: 

1 fy 

lim hd^h = —r^{t) for t > 0 and 2 ; = Zo{t). (1-ld) 

z\zo{t) at 

This condition may be viewed as a Rankine-Hugoniot condition for a viscous shock wave: Since 
(1.1a) has the form of a (nonlinear) continuity equation 

dth -|- dz{Vh) =0 for t > 0 and z > Zo{t), 


(1.2a) 
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where 

V = hdlh (1.2b) 

is the transport velocity of h, by compatibility necessarily (l.ld) holds trne. 

Eqnation (1.1a) is a particular version of the general class of thin-hlm equations 

dth + dz {h"'d^h) =0 for t > 0 and z ^ {h > 0}, (1.3) 

where n is a real parameter. Apparently for n < 0 equation (1.3) has inhnite speed of propagation 
and non-negativity of h is not ensured. On the other hand, for n = 3 (corresponding to the no¬ 
slip condition at the liquid-solid interface) equation (1.3) exhibits unphysical features as well: The 
solution is singular at the free boundary, which does not move unless an infinite amount of energy 
to overcome dissipation is inserted into the system [9,17,26]. Since for n > 3 equation (1.3) is even 
more degenerate, for a model relevant for the motion of fluid Elms (in which contact lines move with 
hnite and in general nonzero velocity) necessarily n G (0,3). The most important cases are n = 1 
and n = 2, corresponding to the lubrication approximation of Darcy’s flow in the Hele-Shaw cell 
{n = 1)^ or the lubrication approximation of the Navier-Stokes equations with a (linear) Navier- 
slip condition {n = 2)^ at the liquid-solid interface, respectively [8,10,22,23,27]. More detailed 
discussions of the literature, also addressing the well-established existence theory of weak solutions, 
can be found in [2,4,15]. 

This note addresses the regularity of solutions at the free boundary and may therefore be considered 
as a contribution to a regularity theory of higher-order degenerate-parabolic equations, an only 
insufficiently explored field. The author hopes that the present study is also relevant regarding an 
existence, uniqueness, and regularity theory for the (Navier-)Stokes equations with a moving contact 
line, an essentially open problem. Here, a thorough understanding of the singular behavior at the 
free boundary can potentially help in the construction of suitable function spaces. Additionally, 
in view of the aforementioned no-slip paradox, the question of how the Navier-slip condition (or 
even general nonlinear slip conditions) de-singularizes the solution at the contact line 
is also of interest from the applied view point. Loosely speaking, the following arguments will 
demonstrate that higher spatial regularity (i.e., regularity in the variable z) goes hand in hand with 
higher regularity in the time variable t. In numerical studies in [28], this observation is used to test 
numerical schemes regarding their precision in resolving the dynamics of moving contact lines. 

1.2. Transformations. We review the transformations discussed in detail in [11, Sec. 1.2, Sec. 1.3, 
Sec. A.l]: Generic solutions to (1.1) are traveling-wave solutions, that is, solutions of the form 
h(t, z) = where x = z — V^t and Vq > 0 is the speed of the traveling wave. By rescaling, 

we may without loss of generality assume Vq = | and one may conclude that in the case of a 
moving infinite cusp a strictly monotone similarity solution is given by^ H'y^^(x) = x^. Considering 
perturbations of this profile, we set 

h{t, Z{t,x)) = x^ ioTt,x>0. (1.4) 

Under the assumption that h(t,z) is strictly monotone in z ioi z > Z^it), equation (1.4) uniquely 
determines the function Z = Z{t,x), thus interchanging dependent and independent variables and 
hxing the boundary to t = 0. The transformation (1.4) is known as the hodograph transform 
(cf. Fig. 1). The related von Mises transform has been applied already in the context of the porous- 
medium equation and the thin-hlm equation with linear mobility in higher dimensions [18,24]. 

^For an extensive well-posedness and regularity theory of (1.3) for n = 1 in the complete wetting case, we refer 
to [13,14,16,18]. The partial wetting case is addressed in [22,23]. 

^We refer to [11] for a well-posedness and partial regularity result of (1.3) in the complete wetting regime with 
n = 2 and to [19] for a result in the partial wetting regime. 

discussion of traveling-wave solutions can be found in [5]. 
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Plugging expression (1.4) into (1.1a) and noting that the chain rule (applied to (1.4)) gives the 
transformations 

dth = —Ztdzh and = Z~^dx^ (1.5) 

equation (1.1) now reads^ 

— ZtZ^^dxX"^ + Z~^dxX^ (^Z~^dxY x^ ~ ^ for t,a; > 0. (1.6) 


Dehning 

and noting that Ft = 


F := Z-^ 

—F'^Zxt, we observe that (1.6) alters to 
xdfF + F‘^xdxX~^ dxX^ FdxFdxFx^ = 0 


for f, a; > 0. 


(1.7) 


By commuting the powers of x with the differential operators d^, we then arrive at the equation 


xdtF + M{F, ■ ■ ■ , F) = 0 foTt,x>0, (1.8) 

where is a 5-hnear form explicitly given by 

Af (Fi, • • • , Fs) := F 1 F 2 F (^D + 0 F 3 (^D - 0 F, (^D + 0 F, (1.9) 

and D := xd^ denotes the scaling-invariant (logarithmic) derivative in space (setting s := Inx we 
have D = dg). We observe Ztw ^ = ^ ^ ~ 1^ -^tw ^ = ^ 1 for fhe traveling-wave prohle (indeed 
A4(l, • • • ,1) ^ = ^ 0), so that by setting 

u:=F-l (1.10) 

in fact utw = 0 and we are lead to study the Cauchy problem 

xdtu + p{D)u = Af{u) foTt,x>0, (1.11a) 

u\t=o = for a; > 0. (1.11b) 

Here we use the following notation: 

• p{() is a fourth-order polynomial 

P(C) := C 2 ^ ~ 4) 2) ^ ~ ^ ^ ^ 2) 2) (f-f 2 ) 


with the irrational root f3 := ^ ^ 0.65. The operator p{D) can be derived from the 

5-hnear form M. (cf. (1.9)) by noting that 

p{D)u = Mil, ■ ■ ■ , !,«) H- V Miu, 1, • • ■ ,1). 

• J\f{u) stands for the nonlinearity of the equation and has the structure 

A/’(u) := piD)u — A4(l -|- u, • • ■ , 1 -|- u). (1.13) 

The structure of (1.11a) is such that the left-hand side is linear in u, whereas the right-hand side 
contains terms that are at least quadratic and at most quintic in {u,Du,--- The form 

of (1.11a) can be guessed from (1.1a) immediately: By (1.4) the linearization of the spatial part 
dz {h^dlh) of (1.1a) scales as x~^ and therefore the spatial operator x~^piD) appears, where p(C) 
has to be of order 4. This implies the space-time scaling x ^ t, unlike a; ~ fL for non-degenerate 
fourth-order parabolic equations. Two of the roots of the polynomial p{() are immediate as well: 
The root — | is directly related to the exponent of the right-hand side of the hodograph transform 

^Here and in what follows, differential operators act on everything to their right, that is, e.g. 

Z-^dxX^ (Z-^dxfxi = Z-^dx (x^ (Z-^dx (z-^dx (z-^dx)))) . 
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(1.4). The root 0 appears as the nonlinear equation (1.1a) has divergence form and this feature is 
preserved in a linearization. In other words, one may understand the occurrence of this root by 
noting that the traveling wave Trw = 1 is simultaneously a solution of the respective linear and 
nonlinear equations, and the equations for u and F yield the same linear operator. In contrast 
to these two roots, the emergence of the roots and —/3 — | in the polynomial p{() in ( 1 . 12 ) 
is a genuinely non-trivial feature, specihc to higher-order degenerate-parabolic equations; In the 
second-order case, p{() would be a second-order polynomial and the only two roots of p{() would 
be immediate by the same arguments as above. An accessible way to understand this uncommon 
phenomenon in the fourth-order case for source-type self-similar solutions using the language of 
dynamical systems theory is explained in [3,12]: There the roots discussed above turn out to be the 
eigenvalues of a linearized dynamical system at a stationary point (representing the contact line) 
and the situation can be reduced due to known regularity results of the corresponding invariant 
manifolds (Hartman-Grobman theorem). 

The rest of the note will be mainly concerned with the analysis of (1.11). 


1.3. Notation. For f,g > 0 we write / <s 9 if there exists a constant G > 0 depending on 
parameters S such that / < Cg. In this case we say that / can be estimated by g or equivalently 
that g bounds/controls /. Furthermore, we write / 9 f 9 and / >5 g. For a non-negative 

quantity r we say that a property is true for r <^s 1 (or r ^>5 1 , respectively) if there exists a 
(sufficiently large) constant G > 0 depending on S such that the property is true for r < C~^ 
(r > G). Then we say that r > 0 has to be sufficiently small (large). If A = 0 or if the dependence 
is specihed elsewhere, we just write f ^ g etc. The space G“((0,cxd)) denotes the space of test 
functions, i.e., the space of all u : (0, cxd) M which are inhnitely differentiable with compact 
support contained in (0, cxd). For a G M, we denote by [aj := max{A; E Z : k < a} the integer part 
(floor) of a. We write |A| for the number of elements (cardinality) of a hnite set A. 


1.4. Weighted L^-norms. For the subsequent results, we introduce a scale of weighted L^-norms 


j. These are given by 


Wa ■ = 




X '^{u{x)y 


X 


e (M(e^)) ds with a G 


(1.14) 


The larger a is, the better the decay of ^(a;) as t \ 0 if |m|^ < cxd. To make this a point-wise 
statement, it is convenient to dehne Sobolev norms 

^ /»00 1 ^ poo 

,« •= {D^u{x)Y — = J2 ('9fM(e"))^cls with A; G No, a G M. (1.15) 

/?—n ^0 ^ n ^—00 


■= 


We also use the notation (•, ■)k,a and (•, •)q, for the induced inner products. Setting t(s) := e~°‘^u (e®), 
it is elementary to see that \u\f^^ ~q, ||T||^fc. 2 (R). In particular \u\^^ < 00 implies u{x) = o (a;“) 
as T \ 0. Note, however, that control of an increasing number of H-derivatives does not lead to 
better regularity properties of u{x) in t = 0 as H is scaling-invariant in x. Essentially the index k 
in the norm ^ controls the interior regularity of m, whereas a yields control on the regularity at 
the boundary a; = 0. The identihcation with the standard Sobolev norms also guarantees that the 
test functions G“((0, cxd)) are dense in the spaces {u locally integrable : \u\f,^ < cxd}. More details 
are contained in [ 11 , Sec. 4]. 


The well-posedness proof of [11] relied on the control of the initial data in the norm 


u 


(0)11|2 _ I (0) 


I k-\-6,—S 


+ 


u 




^+6,(5 


Iq, where 

(1.16) 



6 


MANUEL V. GNANN 


k > 3, {x = 0) (the boundary value of and 5 > 0 is chosen snfficiently small. 

Indeed, the snbseqnent result will use the same norm with snfficiently large /c G N and snfficiently 
small (5 > 0 depending on the valne of Nq. We notice that by qnite elementary arguments (cf. [11, 
Eq. (8.5)]) the Lipschitz norm of the hodograph coordinates Z{0,x) — a; is controlled by this norm, 
i.e., snp^-^o I'M®(x)I < |||w®|||q. Smallness of the Lipschitz norm of Z{0,x) —x ensnres strict 
monotonicity and thus invertibility of transformation (1.4). This was in fact the crucial assumption 
in [18,24]. 

2. The main result 

2.1. A regularity result. For motivating onr main result, we may have a heuristic look at the 
properties of problem (1.11): Since we are interested in a pertnrbative result, that is, and u are 
assnmed to be small in snitable norms (cf. e.g. (1.16)), the precise nnderstanding of the linearized 
problem 

xdtu + p{D)u = f ioT t,x > 0, (2.1a) 

u\t=o = for a; > 0 (2.1b) 

with a right-hand side f(t,x), turns ont the be essential. For a; 1, the term p{D)u will be 
dominant. Considering p{D)u ~ / for a; 1, all solntions of this ordinary differential equation 
(ODE) are given by the sum of a particular solution and a linear combination of solntions to the 
homogeneous equation p{D)u = 0. The solution space of p{D)u = 0 is spanned by a;°, x^, x~^~ 2 , 
and x~ 2 . Clearly, the last two powers are ruled out by compatibility with transformations (1.4) and 
(1.7). Nonetheless, the other two, x^ and x^, generically appear in the expansion of the solntion u 
close to the bonndary a; = 0 as a detailed analysis of the resolvent eqnation shows [11, Sec. 6]. Dne 
to the perturbation xdtU, one may then expect that u = v + x^w, where both v and w are smooth 
fnnctions in x np to the bonndary a; = 0. While it is in principle possible to constrnct solntions of 
this form for the linear problem (2.1), snch a structure is incompatible with the nonlinear eqnation 
(1.11a) that mixes the powers x and x^. What we are instead aiming at is 

u{t, x) = u (t, X, a;^) , (2.2) 

where it = u(t, x, y) is a smooth function in (a;, y) G {(a;', y') G : a;', y' > 0}. Indeed, such a result 
can be expected in view of the analysis of source-type self-similar solutions by the author joint with 
Giacomelli and Otto [12], where a stronger result is shown: There 

h{t,z) =t~^H{x) with. X = zt~^ and H{x) = Cx^ (l + u [x,x^)) , 

where C is a positive constant and u{x,y) is analytic in a neighborhood of {x,y) = (0,0). Not 
surprisingly, as in (1.11) or (2.1), (3 = turns ont to be the root of a polynomial that determines 

the linearized problem for H. 

We denote by 

Kno ■= + I3n2 : (ni,n2) G Nq, rii + f3n2 < iVo} (2.3) 

the set of admissible exponents np to order O [x^°) (cf. Fig. 2). The main result of the note reads 
as follows: 

Theorem 2.1. For any Nq G Nq there exist e > 0, k E N, and 5 > 0 such that if 

< 
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the unique solution u of (1.11) (constructed in [11, Th. 3.1]) fulfills the point-wise expansion 
(cf. Fig. 2) 

u{t,x)= Ui{t)x^ + Ri^fix,t)x^° as a: \ 0 and t > 0, (2.4) 

where the coefficients Ui = ufif) are continuous functions in t > 0, and the correction = 
i? 7 Vo(x, t) is again continuous and for every t > 0 uniformly hounded in x > t). Furthermore, 


and 


Ui{t) = o (t *) as t —)■ oo 


(2.5a) 


Rn^xR) =o{t ^°) as t ^ oo uniformly in x > 0. 


(2.5b) 



Figure 2. Schematic of the admissible exponents i G (cf. (2.3)) in the general¬ 
ized power series (2.4) of the solution m up to i < 4. 

Actually, we are able to prove further regularity properties for the coefficients Ui{t) and the correction 
Rno = RNo{x,t) and also interior-regularity bounds, the presentation of which we postpone to later 
sections (cf. Lemma 3.4 and Theorem 4.1). Note that in [11] it was only proven that u{t,x) = 
uo(t) -\- uy{t)x^ o{x^) as a: \ 0 almost everywhere, where uq = uo(t) is bounded and continuous 

with Uo{t) —)■ 0 as t —)■ cxo and t^~^ui 3 = t^~^up{t) is square integrable with uy{t) = o 
t —)■ oo for a subsequence. 

2.2. Discussion. The result captures the regularizing effect of the degenerate-parabolic equation 
(1.11a). Unlike in the case of the porous-medium equation 

dth — d]h^ = 0 for t > 0 and z e {h > 0}, 
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with a constant m > 1 - the second-order connterpart of (1.11a) for which a comparison principle 
holds and solntions become smooth for positive times [1, 7, 24] - here the solntion only becomes 
smooth in the generalized sense (2.4). 

We emphasize that in the case of partial wetting, that is, withont loss of generality \dzh\ = 1 
a,t z = Zo(t), Kniipfer fonnd that generically snch singnlar expansions do appear as well; The 
solntions turn out to have a polynomial expansion in x and x\n.x for n = 2 [19], and x and x^~^ for 
n G (O, \ {1,2, |, |, [20,21] (cf. also [3] for a discussion of source-type self-similar solutions 

with nonzero dynamic contact angle). In these cases, however, Kniipfer assumes such an expansion 
already for the initial data which also need to fulhll additional compatibility conditions. Hence in 
this case there is no proof available that the solution acquires additional regularity at the contact 
line for positive times, except for a smoothing effect in [20, 21, Cor. 4.3] if one waits sufficiently 
long. This is due to the different techniques in [19, 20] relying on the Mellin transform and a 
suitable subtraction of the singular expansion multiplied with a cut off at a; = cxd: The resulting 
estimates contain terms with different scaling in x and consequently have no distinct scaling in 
time t. Thus it seems impossible to introduce time weights in order to capture the smoothing effect 
of the (non-)linear equation in this setting. On the other hand, it is well-known in the theory of 
non-degenerate parabolic equations (in domains with smooth boundary for instance) that solutions 
typically become smooth for positive times without having to assume well-prepared initial data. 

This (generalized) smoothing effect is the essentially new insight of Theorem 2.1. The fact that we 
need control of an increasing number k + 6 of logarithmic derivatives D the larger Nq, essentially 
amounts to having sufficient interior regularity of the initial data. Since the degeneracy of (1.11a) 
is immaterial away from a; = 0, the interior regularity of problem (1.11) can be treated by standard 
theory. 

We further remark that e has to be chosen dependent on Nq. With our method it appears to be un¬ 
avoidable to assume this dependence^, so that we cannot determine whether the series Ui{t)x^ 

converges. 

We also notice that there are connections of this work with the theory of elliptic boundary problems 
in domains with isolated point singularities at the boundary, for which it is known that singular 
expansions of solutions occur [25]. This is not surprising, as the underlying fluid model, the Stokes 
(or Navier-Stokes) system, has to be solved on a moving inhnite-cusp domain. In the partial wetting 
case, instead, a moving inhnite wedge may be considered and in fact, Kniipfer’s analysis in [19,20] 
strongly relies on this analogy. 

2.3. Transformation into the original set of variables. In [11, Rm. 3.2, Rm. 3.3, App. A] it 

was explained that analogous leading-order expansions also hold for the function h and the velocity 
V (the vertically-averaged horizontal velocity within lubrication theory, cf. (1.2)). Indeed, in the 
expression V = hd^h for the velocity (cf. (1.2b)), we may employ the hodograph transform (1.4), 
i.e., h = x^ as well as dz = Fdx (cf. (1.5)&(1.7)), so that 

V = ^x^FdxFdxFx^ = M(F, F, F) (2.6a) 

with 

M{F,,F,,Fs) := (^D - 0 F^ (^D + 0 F^. (2.6b) 

^As we will observe in Section 3, the coercivity constant of p{D — Nq) with respect to the inner product (•, da 
vanishes as a A^o or a \ A^o — !> and as Nq ^ oo the admissible exponents in the interval {Nq — 1, A^o) become 
dense. Thus the constant in the maximal-regularity estimate (3.41) blows up as Nq oo. 
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Equations (2.6) can be used to derive 

V = A{t) + B{t,x) + C{t,x) (2.7a) 

with 

7l = -^(l + Mo)^ (2.7b) 

S = 3(1 + uofMsym{u - Uo, 1,1) = ^(1 + uofp{D){l + uq) where p(C) 

C = 3(1 + Uo)Msym{u -Uo,U- Uq, 1) + Msym{u -Uo,U- Uo,U- Mq)- 
With help of expansion (2.4), equations (2.7) upgrade to 

V(t, Z(t,x)) =--{1 + Uo(t))^ + Vi(t)x^ + O [x^°) as T \ 0 and t > 0, (2.8) 

iSAjVg 

l<'i<7Vo 

where the coefficients Vi(t) and the correction O fulhll decay estimates as in (2.5). From (2.8) 

we can also read off the velocity Vo(t) and the position Zoit) of the contact line as 

l/o(t) = ^ (1 + Mo(t))^ and Zo{t) = zo-^[ {1 + uo {t')f dt', 

8 8 Jo 

where ; 2 o ^ ® is a free parameter. 


— ( C + /^ + A ) iC — l 3 ), 

(2.7c) 
(2.7d) 


Finally, we can derive an expression for the expansion of the him height h and the velocity V in 
the vicinity of the contact line in terms of the original variables t and Employing (2.4) in (1.7) 
in conjnnction with (1.10), we conclude 

Z^ = {1 + uo{t))~^ + {Zx)i{t)x^ + O [x^°) as T \ 0 and f > 0, 

iSAjVg 

/3<i<No 


where the coefficients and remainder fnlhll estimates analogous to (2.5). Integration of this expan¬ 
sion gives 


X 


Z{t,x) - Zo{t) 

( 


X 


(l-hMo(t)) ^ 


V 


\ 


(1 -|- i) ^ {Z3;)i{t)x^ + O 

/3<i<No 


{X 


No 


) 


/ 


as T \ 0 and t > 0. (2.9) 


Inversion yields 

T = (1 Uo{t)) X 


/ \ 

1 + ^ Cj(t)7* -f O {x^°) 

ieAjvo 

\ l3<i<No / 


as 7 \ 0 and t > 0, 


( 2 , 10 ) 


with coefficients Ci{t) and remainder O [x^°) obeying estimates analogous to those in (2.5). Utilizing 
expansion (2.10) in the hodograph transform (1.4), we end np with 


h(t, z) = x^ 1 + 


Ui{t)3P + O [x^°) I 

iGKjy / 


as 7 \ 0 and f > 0, 
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where the Ui(t) are (at least) continuous in time t, x := z — Zglt) (cf. (2.9)), and decay estimates 
as in (2.5) for the Ui(t) and the remainder O (x^°) hold true. We leave it up to the reader to 
directly relate the coefficients Ui(t) to the Ui(t) appearing in (2.4). Furthermore, we can also derive 
an expansion for the velocity V = V{t,z) by using (2.10) in (2.8): 

V(t, z) = —- (1 + uo(t))^ + Vi(t)x^ + O (5;^°) as 5; \ 0 and t > 0, 

1<2<A^0 

where again x := z — Z^it) and the Vi(t) and the remainder term O (5;^°) meet decay estimates 
as in (2.5). We remark that it is open, whether a similar expansion for the velocity occurs for the 
(Navier-)Stokes equations on a moving inhnite-cusp domain. 

2.4. Outline. The paper consists of two parts: the linear theory (discussed in Section 3) and the 
nonlinear theory (contained in Section 4). After recalling some of the basic notions on coercivity 
and maximal regularity in Section 3.1, Section 3.2 deals with the formal structure of the linear 
degenerate-parabolic equation (2.1). Here we demonstrate that by applying an appropriate combi¬ 
nation of scaling-invariant operators H — 7 with 7 G M to the linear equation ( 2 . 1 a), we are able to 
derive maximal-regularity estimates for ( 2 . 1 ) that control the singular expansion of m at a; = 0 to 
arbitrary orders (cf. Section 3.3). Since all higher-order equations for u are scaling-invariant, the 
corresponding maximal-regularity estimates have a distinct scaling in x. Thus, when multiplied with 
an appropriate time weight, we can combine them to a quasi scaling-invariant maximal-regularity 
estimate for u in terms of the initial data in the guasi-minimal norm (1.16) and the right-hand 
side /. These arguments are made rigorous in Sections 3.4 and 3.5 (cf. Proposition 3.7) without 
going into all details. In particular Section 3.5 is not essential for the understanding of the main 
ideas. 

Finally, in Section 4 we prove our main regularity result. Theorem 4.1, from which Theorem 2.1 
follows as a special case. The proof strategy is standard and requires two ingredients: maximal- 
regularity estimates for the linearized problem ( 2 . 1 ) and the factorization of the nonlinearity J\f{u) 
(cf. (1.13)) given in Proposition 4.2. The latter is in fact the non-trivial part of the proof (cf. Sec¬ 
tion 4.2) and requires detailed estimates that rely on the symmetry properties of the multi-linear 
form Ai (cf. (1.9)). 


3. The linear problem 

3.1. Coercivity and parabolic maximal regularity. Again, we briefly repeat some of the linear 
theory in [11]. Consider the linear problem 

xdtu P{D)u = f ior t,x > 0, (3.1a) 

u\t=o = (3.1b) 

which is structurally the same as ( 2 . 1 ) but with a general fourth-order polynomial P{C)- We assume 
that the zeros 71 < • ■ • < 74 of P{C) are real. Then we know from [ 11 , Prop. 5.3] that there is 
a range of weights a - which we may call coercivity range - such that P{D) is formally coercive, 
i.e., {u, P{D)u)a Wlla ^ ^ Sufficient criterion (which can be elementarily 

computed) is 

2 / \ 

a e (-cx),7i) n (72,73) n (74, oo) and {a - m{'y)f < , (3.2) 

where 771 ( 7 ) •= i (piean of the zeros 7 ^) and := j X]j=i( 7 i ~ (variance of the 

zeros 7 j). In the particular case of P{D) = p{D), (3.2) yields that the coercivity range contains the 
interval (— 1 , 0 ). 
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Now suppose that a G M is in the coercivity range of P{D). Then at least formally by quite 
elementary arguments (cf. [11, Sec. 2 , Sec. 7.1]), we can derive a differential version of a maximal- 
regularity estimate for (3.1a) that reads 

^ |w|?+ 2 ,a-i + l/l^ > ^here ^ e No. (3.3) 

Indeed, in (3.3) derivatives /, • • • , D^f control u, • • ■ , in the same norm, which is the maximal 

control in space one can expect as (3.1a) is fourth-order in D. The additional control of the time 
derivative dtu, • • • , D^dtu (with reduced weight due to the degeneracy in (3.1)) can be obtained by 
using control of u, • • ■ , and the fact that (3.1a) is fulhlled. This also yields control of the 

trace ^ |'u |£_,_2 q,_i by interpolation. We refer to Section 3.5 for more details. 

By multiplying (3.3) with a time weight P (where a > 0), we obtain the integrated version of (3.3), 

i.e.. 


sup t"" \ufi+2,a-l + + Iu|?+4,a) 

oOO POO 

Jit,O Ao|“*°’|,+2,0-1 + / + l“lm,a-l where No. 


(3.4) 


The purpose of introducing time weights is two-fold: On the one hand they will enable us to 
prove the decay estimates as f —)■ cx) in (2.5) for the coefficients u* and the remainder While 
this would be irrelevant on a hnite time interval, they also make it possible to prove regularity 
immediately after time t = 0, whereas without them the time after which regularity is obtained is 
unknown (cf. [20,21, Cor. 4.3] for a similar case). 


Estimate (3.4) will be the basis of all linear estimates in the sequel. 


3.2. The formal structure of the linear equation. As pointed out in [ 11 , Sec. 2 ], just applying 
maximal regularity of the form (3.4) with a = 0 to the linear equation (2.1a) is not sufficient in 
order to obtain well-posedness of the corresponding nonlinear problem (1.11). This is so, since only 
negative weights a are admissible (viz. in the coercivity range of p{D)) and hence no control of the 
boundary value uo or the sup-norm sup^ \u{t, x) \ can be achieved. On the other hand, as products 
of up to hve factors in {u,Du,--- ,D^u} appear in the nonlinearity Miu) (cf. (1.9) and (1.13)), 
the control of sup^a-^g \u{t,x)\ appears to be necessary for proving well-posedness by a contraction 
argument. In order to circumvent this problem, it was convenient to apply p{D — 1) to the linear 
equation (2.1a) that, by using the commutation relation Dx = x{D + 1), transforms into 

xdtv^^'^ + p{D — l)u*^^^ = for t,x > 0, (3.5) 

where := p{D)u and := p{D — 1)/. Since the coercivity range has translated to the positive 
interval (0,1) (cf. e.g. (3.2)), one obtains better control on the regularity of v at the boundary 
T = 0. This is not surprising as in view of (2.4) we expect u{t, x) = u^it) + Uj 3 {t)x^ -f 0(a;) as t \ 0 
and t > 0 and the powers x^ and x^ are in the kernel of p{D), whence v^P{t,x) = 0{x) as t \ 0 
and t > 0. Furthermore, by compatibility f(t,x) = 0(x) as x \ 0 and t > 0 and therefore also 
gP{t,x) = 0{x) as T \ 0 and t > 0 . As a second step one can then retrieve regularity information 
on u from regularity information on using elliptic estimates for the operator p{D) (cf. [11, Sec. 2, 
Lem. 7.2]). 

Before reviewing the arguments, we point out the limitations of the ansatz: As a natural second 
step, we apply the operator p{D — 2) to equation (3.5) and obtain 

xdtv^'^^ + p{D — = g^"^^ for t,x > 0, 


(3.6) 
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with := p{D — and := p{D — 2)g^^\ Following the argumentation above and not¬ 

ing that the coercivity range of the operator p{D — 2) contains the interval (1,2), we apply the 
maximal-regularity estimate (3.4) to (3.6) and seemingly obtain even better control on the bound¬ 
ary regularity of which formally suggests v{t,x) = O as a; \ 0 and t > 0, where e > 0 

is arbitrarily small. Apparently, such a claim is too strong as generically also the term x"^^ appears 
in the expansion of u (cf. Section 2.1), this term is not in the kernel of p{D — l)p{D), and thus in 
general n(f, x) = O as a; \ 0 and t > 0. 


In order to work around this problem, we set 

h := {ni + /3n2 : (ni, 112 ) e Nq, 1 < tii + /3n2 < 2} \ {1 -F /?} = {2/3, 3/3} 

and apply the operator = (-D —2/3)(Zl —3/3) to (3.6). Setting := (11*6/2~ 

and := (11*6/2~ *)) (3.6) turns into 


xdt'w^'^^ + p{D — 2)w^'^'^ = -|- xq 2 {D)dtV^‘^^ for f, a; > 0, (3.7) 

where q 2 {D) = 11*6/2“ 11*6/2~ —2D + 5/3 — 1 is a polynomial of degree I/ 2 I — 1 = 1. 
The additional term xq 2 {D)dtV^‘^^ appears as the commutator of 11*6/2multiplication 
with X does not vanish. Nevertheless, (3.7) is structurally advantageous compared to (3.6) as now 
we may indeed expect w^‘^\t,x) = 0{x‘^), r^‘^\t,x) = 0{x‘^), and xq 2 {D)dtV^‘^\t, x) = o{x‘^) as a; \ 0 
and f > 0. Hence applying the maximal-regularity estimate (3.4) with a G (1,2), we obtain better 
control on (assuming a > 0): 


t>o 


/+2,a-i 


1 + 


x2cr 




£,a—l 


+ \w 


( 2)1 


I £+4,0; 


dt 


^£,a 


f-\r^^'>\l^dt + 2a +-1 


w 


( 2)1 


l/+2,o-i 


idt+/ \q2iD)d,v^X,a-i^^- 


(3.8) 


Yet, the additional term /g°°+ \q 2 {D)dtV^‘^^\^^ ^_^dt in (3.8) has to be treated: Since a — 1 G (0,1) 
is in the coercivity range of p{D — 1), we can apply maximal regularity of the form (3.4) to the 
time-differentiated version of (3.5), i.e.. 


a;( 9 +d) — l)dtv^^^ = dtg^^^ for f, a; > 0, 


(3.9) 


which leaves us with 

sup+ |i9tnd)| 


i>0 


< 


e'+2,a- 


3 + / t 


^ £', 0—2 




£'+4,0-1 


dt 


£',o / t'" dt + 2a +-1 i^^^(i) 


(3.10) 


£'+ 2 ,o-§ 


Assuming £' > £ -|- I/ 2 I — 5 = £ — 3 and noting that trivially \q 2 {D)dtV^‘^^\‘^^ l £'+4 o’ 
combination of (3.8) and (3.10) yields 


sup t 
t>o 


2(7 


w 


( 2)1 


|£+2,q- 


1 + 


£'+ 2 ,Q-j 


+ (la.u.<='lh_. + + \ayX..-. + |a.-'Oh4,„-0 

«',a [ «"“(k'"’lh+|9i9“'l/,„_i)d« 


dt 


'0 

< 


(3.11) 


2a / t 

Jo 


2(7-1 


w 


( 2)1 


I£+2,q-4 


1 + dt) dt. 
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Now the solutions and are estimated in sufficiently strong norms by and in 

respective norms and the integral 

+ (3.12) 

The norms in (3.12) have reduced spatial and temporal weights and thus it is possible to absorb 
these terms by lower-order estimates®. We will detail the arguments in Section 3.3. Notably it 
seems unavoidable to combine higher-regularity estimates in space with higher-regularity estimates 
in time as opposed to the case of the thin-hlm equation (1.3) with n = 1, that is, the lubrication 
approximation of the Hele-Shaw cell [14, Sec. 8, Sec. 9]. This has already been observed by Kniipfer 
in [19,20] for the partial wetting case using rather different techniques. 

Before addressing the issue of dealing with the lower-order terms in (3.11), we will hrst systematize 
our observations: In order to obtain better control on the solution, we apply the operator p{D — 3) 
to equation (3.7). Thus we arrive at 

xdtv^^^ + p{D — 3)n*^^^ = + xq 2 {D)p{D — 2)dtV^‘^^ for t,x > 0, (3.13) 

where we have set := p{D — and g^^^ := p{D — 3)r^^b Again, we cannot expect to have 

x) = 0{x^) as T \ 0 and t > 0 as the set 

-^3 + /^'^2 ( (RI) R 2 ) ^ Nq; 2 < Til -\- /3n2 < 3| \ {2 -f- /?} 

is non-empty (cf. Fig. 2). Applying nie/ 3 (-^ ~ *) (3-7), we obtain 

xdtw^^^ + p{D — 

= -f xq‘i{D)dtV^^"’ + xq2{D)p{D - 2) | ]^(T) - i -h 1) 

V*e/3 

where again ;= ([ 1 ^ 6 / 3(71 - i)) r^^'> := ( 11 * 6 / 3 “ *)) ^3(A>) = 11 * 6 / 3 - 0 “ 

— / -f 1) is a polynomial in D of degree I/ 3 I — 1 that originates from the commutator of 
X and 11 * 6 / 3(77 ~ *)• Then we make the rather trivial observation /2 C Is — 1 := {/ — 1 : / G Is} 
(cf. Fig. 2). By exploiting 


for f, a; > 0, 


(3.14) 


p{D-2)ll[{D-t + l)]d,v^^^=[ J] (D-z) p{D-2)ll[iD 


t\v 


( 2 ) 


\ieI 3 


,*e(/3-i)\/2 


vie/2 


=*h(2) 


n (22 - >) 

,*e(/3-l)\/2 


.(3) 


and setting 

qs{D) := UD) + q 2 {D) J] (77-*), 

*e(/3-i)\/2 

which is a polynomial of degree less or equal to |Is| — 1, we can rewrite the last two terms in (3.14) 


and obtain: 


xdtw^^^ + p{D — -\- xq 3 {D)dtV^^^ for f, a; > 0. 


(3.15) 


6 


Here “lower order” is meant in the sense of using norms with lowered weights and the same number of time 


derivatives. 
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As (3.15) is structurally the same as (3.7), this procedure can be iterated, and we arrive at the 
following set of equations 

{xdt + p{D - n))d^w^^^ = + xqn{D)dr^h^^'> for f, a; > 0, (3.16) 

where we dehne the sets of indices 


In ■= {^1 + /d ^2 : {ni,n 2 ) e Nq, n - 1 < Ui + /3n2 < n} \ {n - 1 + /9}, (3.17a) 

Jn ■= {ni + /3n2 : (ni, 742 ) G Nq, 0 < ni + /3n2 < n} \ (N U (No + /3)) = (3.17b) 

introduce the functions 


V 


(n) 


W 


(n) 


'n-1 


JJ p{D - n') 

rL'=0 

n—1 

JJ p{D - n') 


\n'=0 



(3.18a) 

(3.18b) 


ri^) ■ — 


np(^-n') n(^-*) 


\n^=l / \i£Jn / 

and denote by qn{D) a polynomial of degree less or equal to |/„| — 1. 
m G No are arbitrary. 


(3.18c) 

The numbers n G N and 


3.3. Heuristics for parabolic maximal regularity. In this section we systematize the ideas of 
the previous section, leading to estimate (3.11). Throughout the section, all estimates may depend 
on N, n, m, a, or S. Applying the maximal-regularity estimate (3.4) to equation (3.15), we obtain 


sup 

t>o 


k(n,m,a')+2,a'+n—^ 


+ [ j2{a+n+m)-Z 


L o 9 + 

K(n,m,Q:'),Q:'+n—2 ^ 


fc(n,m,Q:')+4,Q:'+n—1 


dt 


< 5- 


2(a+n+m—l),l 


am.,/”) 


k{n,m,a')-\-2,a'-\-n—-. 


+ / (3.19) 


+ (l-(5„,l) / t2(a+n+m)-3 g^+ly{ 


n)\ 


I /c(n,m,Q;^) + |7n| —2 


dt 


+ (2(a + n + m) - 3) / 


k{n,m,a')+2,a' +n—^ 


3 dt. 


where we assume and use the following: 


• We take a hnite number of weights a G [0,1] and we choose 5 > 0 sufficiently small such 
that 

A a' := a ± 5 G (0,1) if a G (0,1), 

A 5 G (0,1) and thus also 1 — <5 G (0,1). 

Thus a' is in the coercivity range of p{D — 1) (i.e., o; -|- n — 1 G (n — 1, n) is in the coercivity 
range of p{D — n)). 

• (5 > 0 has to be chosen small enough such that if Oi < 0 : 2 , then also ai+6 < a 2 — S. Further 
smallness conditions on S will be specihed when necessary. 

• We need to assume 2 (q; -|- n -|- m) — 3 > 0 so that all time weights are integrable at t = 0 
(the last term in (3.19) vanishes for 2 (q; -|- n -|- m) — 3 = 0). 

• The indices k{n, m, a) G Nq, determining the number of H-derivatives in the norms appear¬ 
ing in (3.19), will be chosen later. 
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Figure 3. Schematic: absorption of remnant terms I. Each node + corresponds to an 
estimate of the form (3.19). The displayed arrows visualize the absorption mechanism 
for = 4, that is, the remnant at the base of the arrow (forming the second but 
last line in (3.19)) is absorbed by the corresponding estimate (3.19) at the tip of the 
arrow under the assumption that (3.22) holds true. 


The specihc choice of the weights a will be explained further below. Choosing proper weights turns 
out to be essential for obtaining control on the coefficients Ui{t) of the generalized power series (2.4) 
of u and for being able to absorb the last line of (3.19). 

Note that estimate (3.19) itself is insufficient as the solution u still appears on the estimate’s right- 
hand side in the last two lines. 


3.3.1. Absorption of remnant terms I. We will hrst concentrate on absorbing the second but last 
line of (3.19) by exploiting the additional time regularity: Therefore we start with estimate (3.19) 
with n = N and m = 0. For > 1, the term 



^2(a+7V)-3 |^^^(iV)| 


2 

fc(Ar,0,a')+|/iv|-l,a'+Y-2 


dt 


has to be absorbed. This term can be estimated by the left-hand side of (3.19) for n = N — 1 and 
m = 1, provided that the indices k{n,m,a') obey k{N — 1,1,0;') > k{N,0,a') + |/Ar| — 1. Then 
indeed 




k{N,0,a')+\lN\-l,a'+N-2 


< 

r\j 




k(N-\-N—2 


Next, supposed that N > 2, the term 



.2{a+N)-3 |a2 (A-1)|2 

I i lfc(A-l,l,a')+pjv_i|-l,a'+Y-3 


dt 


has to be absorbed, which can be achieved by combining it with estimate (3.19) for n = N — 2 and 
m = 2, provided that k{N — 2, 2, o') > k{N — 1,1, o') -|- |/Ar-i| — 1. Apparently, this procedure can 
be iterated (cf. Fig. 3) and by merely summing (3.19) for indices 


(n,m)e{(iV,0),(iV-l,l),...,(l,iV-l)}, 
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we end up with 


7V-1 




t>0 


m=0 

N-1 


k{N—m,m,a')+2,a'+N—m—i 


poo ^ 

^ / ^2{a+N)-3 J2 \d^+^ 

Jo _ _r, 


W 


(A-m) |2 


m=0 

N-1 


^ / ^2ia+N)-3 


m=0 


I k{N—m,m,a') ,a'+N—m—2 

2 

k(N—m,'m,a')-\-4^a' -\-N—m—l 

/>rv-> —1 


dt 


dt 


(3.20) 




W 


( 1 ) 

|t=0 


/ ^2(a+A)-3 E 


lA3(l,0,i=b5)+2,=b5 Jq 

poo 

+ {2a + 2N-3) / t2(«+^)-4 ^ 


k{N—m,m,a'),a'+N—m—l 


dt 


m=0 


I k{N—m^m^Q')+2,a'+N—m— ^ 


3 dt, 


m=0 


provided that we have 

k{N — m — l,m + l,a') > k{N — m, m, a') + \lN-m\ — 1 for m = 0, • • • , iV — 2 

and 

2a + 2iV-3>0. (3.21) 

For later purpose we require the stronger assumption 

k{N — m — 1, m + 1, a') > k{N — m, m, a') + |/A-m| for m = 0, • • • ,N — 2. (3.22) 


3.3.2. Absorption of remnant terms II. Estimate (3.20) is still insufficient as yet the solution appears 
on the right-hand side of the estimate (forming the last line). The absorption of this last line is 
indeed more complicated and demands to specihcally choose the weights a. In order to understand 
the choices made, we first make additional comments on elliptic maximal regularity: In fact, one 
can get control on u from control on the functions 

Proposition 3.1. Suppose that k G No, G M \ K^o, and u : (0, cxd) —)■ R is smooth satisfying 

D^u{x) = Uii^x^ + o{x^) as X \ 0 (3.23) 

i&Kg 

for all i = 0, ■ ■ ■ , k + \K^\ (ef. (2.3) for the definition of Kg). Then for any polynomial Q(() = 
n^i(C ~ Ce) with real zeros Cu •' ‘ Xm sueh that Kg C. {Cu • ‘ ‘ > Cm} and g ^ (Ci, • • • , Cm}, we have 

<k,XQ{D)u\,g. (3.24) 

A proof in a similar case can be found in [11, Lem. 7.2]. It is an immediate consequence of a version 
of Hardy’s inequality: 

Lemma 3.2. For any w G C°°((0, cx))), 7 , G R with 7 7 ^ ^?, < 00 , and w{x) = o{x^) as 

x\0 (x X 00 ) if 'y < g X > g), we have 


u 


-E 

i&Kn 


UiX 
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For a particular admissible exponent i G K^o \ {0}, we may choose the weights i ± 6. Both i + 6 
and i — 6 appear in the coercivity range of p{D — [ij) or p{D — [ij — 1) (these are the intervals 
(W “1) W) (W) W + !))• Having control on 


t 


2i-l 


w 


(Ld)l 


I k{ [ij ,0,i— [ij +1—5)+4,^—(5 


dt if i G No 


or 


we obtain control on 


^2i+l |y^(Ld+i) 


k(\_i\ +1,0,2— [2j±5)+4,2±(5 


dt else, 


poo 

L 

U — UjX^ 


j&Ki 


dt for t G No 


or else 


and 


poo 

L 

U — UjX^ 


j&Ki 


k{ [ij ,0,2— [ij +1—(5)+4+4 [ij +1 | ,2—5 

2 


dt, 


poo 

L 

u — UjX^ — nix'" 


j&Ki 


[ij +1,0,?— [?J —(l)+4+4( [?J +1)+ 1 J I ,?—S 
2 


dt, 

^ fc([ij +1,0,?— [?J +5)+4+4( Ld Tl) +1 1 ,?+^ 

respectively, by using elliptic maximal regularity given by (3.24). It is quite apparent that, by 
applying the triangle inequality, the last three terms yield control on the coefficient Ui{t) of the 
form: |l^?(t)|^dt. As we have already noted in the introduction, for a = | and N = 1 

the first term in estimate (3.20) also yields control on supj>o l'*^o(OI^ analogous estimates of 
higher-order coefficients are possible. We will postpone the details on how to extract further control 
on the coefficients to Section 3.4 (cf. Lemma 3.4)^. For the moment we just note that for controlling 
the singular expansion of u, it is convenient to use spatial weights i 6 with i G K^o \ {0} in our 
norm. 


Without further ado, we start with the choice of weights by recalling that 

In -1 C/at — l = {i — 1: iG In} (cf. Fig. 2). (3.26) 

We £x Aq ^ hi and are aiming at controlling expansion (2.4) up to order O In view of (3.26) 

it is reasonable to view the set of indices ^ as a subset of [0,1] x {1, • • • , Ad}. For each (a. A) G 4, 
we may use estimate (3.20) where a' = a±S G (0,1). We distinguish between two classes of weights: 

(a) By the above considerations on estimating the coefficients and in view of the fact that (3.26) 
holds true, we start by including (a. A) with 

A Q; G (hjVo — A}) -|- 1) U {/3} and A = 2, • • • , Ag, 

A a G {{Ino - Ag -d 1) U {/S}) n (|, l) and A = 1, 

A as well as (a. A) = (0, A) with A = 2, • • • , Ag 
A and {a, A) = (1, A) with A = 1, • • • , Ag — 1, 
in the set A. Thus we are already able to control all coefficients Ui{t) with Nq > i > I. 

(b) In all these cases, we need to be able to absorb the respective remnant terms forming the 
last line of (3.20). Since the weight in the norm is shifted by this requires to include 
{a, A) with 

A a G ((+0 - Ag -d I) U {/3 - |}) n (O, I) and A = 2, • • • , Ag, 

^Indeed, sufficiently strong estimates on them are essential for proving appropriate estimates for the nonlinearity 
(cf. Proposition 4.2). 
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A a e {Ino - ^0 + I) n (|, l) and iV = 1, • • • , iVo - 1, 

A as well as (a, N) = iV) with N = 1, ■■■ , Nq, 

in A. 

Now that we have chosen the set of weights A, we need to ensnre that the absorption mechanism for 
the last line of (3.20) works. This requires some additional conditions on the number of derivatives 
k{n,m,a'). We note that since (3 is irrational, by choosing <5 > 0 sufficiently small we have a' ^ ^ 
for all (a, N) G A. For (a, N) ^ A such that a + iV > 2, we need to distinguish between three cases 
(cf. Fig. 4): 

(a) If a' > we can absorb the remnant in the last line of (3.20) through 


j.2(Q+Af)-4 ffn AN-m) ^ 

t k{N—m,m,a')+2,a'+N—m—A 


3 dt 


< 

r\j 


.2(a-i+A)-3 I om (A-m)|2 i, 

I 1 \ k(^N—m,m,a' — A^+4,a' — ^+N—m—l ’ 


where the second line of the inequality appears on the left-hand side of (3.20) with a replaced 
by a — i. This requires that the indices obey 

k — m, m, a' — >k{N — m, m, a') — 2 if a' G . (3.27) 

Indeed one may verify that by construction in all such cases (a — |,iV) G A. 

(b) If a' < I and N — m > 2, we can absorb the remnant term in (3.20) through 


j.2(o-|-A)-4 f)m JN-m) ^ 

t k{N—m,m,a')+2,a'+N—m— 


3 dt 


/ CXD 

,2((Q+|)+(Ar-l))-3 I om ((A-l)-m) |2 i, 

1 1 I A:((A—1)—m,'m,o' + i)-|-4,o'-|-5 + (W—1)—m—1 ’ 

where the last line of the estimate appears on the left-hand side of (3.20) with a replaced 
by a -|- I and N replaced by iV — 1. In view of (3.18), this requires the constraint 


k — m — 1, m, a' + -^ > k {N — m, m, a') -|- 2 -|- |/ 7 v-m| if cn' G ^0, . 


(3.28) 


Again, by construction we have (a iV — l) E A. 

(c) If a' < I and iV — m = 1, necessarily m > 1 and we can estimate the remnant in (3.20) by 


^2(a+A)-4 2 


< 


/*oo 

/ ^2((a+|)+A—l)—3 

dim-D+lwd) 

Jo 



A:(l,m—l,a'+i),Q'+| —1 


dt. 


Here, the second line in the estimate is controlled in (3.20) with a + ^ instead of a and m 
replaced by m — 1. Yet, the absorption only works if the indices obey 


k [ l,m — l,a' + > k {1, m, a') + 2 if a' G ^0, - 


(3.29) 


The above argumentation shows that we can restrict our considerations to the cases in which 
q; G (^, l), iV = 1, and m = 0 since in the case a = |, the remnant term in (3.20) disappears. The 
remaining terms can be treated by applying an anisotropic version of Hardy’s inequality: 
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n = N — m 
4 

3 

2 

1 


+ + 

(g) a 


|b 

+ + 




+ + 


{gx-(g) 

|b 
c 1 
<— (g) a 



Figure 4. Schematic: absorption of remnant terms II. Each node + corresponds to 
an estimate of the form (3.20) (with n = N — m). The displayed arrows visnalize 
the absorption mechanism starting with N = 7, m = 3, and a' > |. The symbol 
(8) denotes a shift of a by (keeping n and m hxed). The remnant at the base of 
the arrow (forming the last line in (3.20)) is absorbed by the corresponding estimate 
(3.20) at the tip of the arrow under the assumption that (3.27), (3.28), and (3.29), 
respectively, holds true. 


Lemma 3.3. Suppose v : (0, oo)^ —)■ M is smooth, i G Nq, and a G (|, l). Then 



\(.,a-\±S 


df < 








(3.30) 


A proof can be found in [11, Lem. 7.5]. The remnant in (3.20) for {N,m) = (1,0) is of the form 




2 


dt 


and can be absorbed by estimate (3.20) with a = | and iV = 1, that is. 



2 

fc(l,0,a')+2,a'- 






2 

k,-^±5 


dt + 



2 

/c+4,^di5 


dt, 


(3.31) 


where we write k := k (l,0, | ± (assuming that the values for + and — coincide), provided we 
have 


k > A: (1,0, a') + 2 



and a 

2 


(3.32) 


By summing over all estimates (3.20) with {a,N) G A and a' = a ± 6 E (0,1), and fulhlling 
conditions (3.22), (3.27), (3.28), (3.29), and (3.32), we obtain 


\u\ 


< 

rsj 


\u 


( 0)1 


I *,0 



(3.33) 




2 

*,0 



2 

fc+2,i-<5 




2 

fc+2,^+(5 


where 


with :=p{D)u^°'> 


(3.34) 
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is the norm for the initial data, 


Af-l 


Ml := 


5 ; 

- t>0 ±-l 


12 


(a,N)eA 
Q:^=adiiG(0,l) 


m=0 


\k(N —m,m,Q;')+2,Q:'+A^—m—^ 


POO N—1 

/ ^2(a+A)-3 Y 

■\T\^ A 0 ^_n 


(A-m) |2 


1 k(N—m,m,a'),a'-j-N—m—2 


(a,N)GA 

a^=aiLSG(0,l) 


m=0 

N-1 


/oo ^ 

,iv ifc^ m =0 




lk(N —m,m,a')+4,a'+N—m—l 


dt 


dt 


(a,N)eA 

a'=aikS£(0,l) 


is the norm for the solution u (with dehned in (3.18b)), and 


2 


POO ^ 1 

^ ^ / ^2(a+7V)-3 ^ ^ l^m^(iV-m) 

,N)&A m=0 


fc(A —m,m,a'),Q;'+A^—m—1 


dt 


(ck,7V)G^ 
Q:^=adi(5G(0,l) 


is the norm for the right-hand side / (with dehned in (3.18c)). 
Now we make a further assumption, that is, we assume 


k{n,m,a') is constant for a' G ^^’ 2 ^ ’ ^ (^ 2 ’^) respectively. 


except for half integers a E {0, -,1 


(3.35) 


(3.36) 


(3.37) 


One may verify that thus still conditions (3.22), (3.27), (3.28), (3.29), and (3.32) can be satished 
(see below). 


3.3.3. Maximal regularity for the linear equation. Further applying elliptic maximal regularity given 


by Proposition 3.1, we infer that the norms HI-IHq and 
respectively, are equivalent, where HI-IHq is given by (1.16), i.e.. 


I*, 0 ’ 


and IIHII^, as well as HI-IHj^ and 




\u 


Will := mW 


I /c+6,—5 


+ 




fc+6,5 


the norm lll-IH for the solution u is given by 


N-l 


M : = 


t>o 


(a,N)&A 

q;'=o±5s(0,1) 


m=0 


dTu - 


Y. 


N-l 


Y I i2(„+K,-3 


{a,N)&A ' 

a'=Q±5s(0,l) 


m=0 


N-l 


i<a'+N—m—^ 


ar^u - Y1 


d'^Ui 

"df^ 


-X 


i{N—m,m,a')+2,a'+N—m—^ 


^m +1 


Uj 


i<a'-\-N—m—2 


d^m+l 


■X 


dt 


((N—m,m,a') ,a'+N—m—2 


/ oo ^ 

^2(a+A)-3 Y 

,iv )fc.A rn=0 


(a,N)eA 
a^=Q:diiG(0,l) 






i<a' -\-N—m—1 


d^Ui , 

dt^ ^ 


dt 




(3.38) 
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with (cf. (3.17) and (3.18)) 


£(n, m, a') := /c(n, m, a') + \ Jn\ + 4n, 


(3.39) 


and the norm lll-IH^ for the right-hand side / reads 


_ noo N—1 

^ / ^2(a+7V)-3 ^ 

,N)&A m=0 


(a,Y)eyl 
a^=Q:diiE(0,l) 


sr 


/- E 

m —1 




-T 


dt. (3.40) 


f(A—m,m,o'),Q;'-|-A—m—1 


Consequently, (3.33) turns into the maximal-regularity estimate 


u 


< 


u 


( 0)1 



(3.41) 


For convenience, we summarize the conditions on the numbers a'): i{n,m,a') is constant 

for all a' G (O, and a' E (|, l), respectively, except for a E {O, |, l}. Furthermore, the following 
inequalities (through (3.39) equivalent to (3.22), (3.27), (3.28), (3.29), and (3.32)) must hold: 


i{N — m — l,m + l,a') > i{N — m, m, a') — 4 ior N — m > 2, 
i — m, m, a' — > i{N — m, m, a') — 2 if a' G 1^ and a 4- iV > 2, 

£ — m — 1, m, a' + -^ > £ {N — m,m, a') — 2 if a' G ^0, , N — m > 2, 

£ ^1, m — 1, a' > £ (1, m, a') 4-2 if a' G ^0, and m > 1, 


'1 \ 1 

k > £{1, 0, a') — 2 if a' G ( -, 1 ] and a ^ 


(3.42a) 

(3.42b) 

(3.42c) 

(3.42d) 

(3.42e) 


It is apparent that conditions (3.42) can be fulhlled and non-negativity of k (n, m, a') can be ensured 
(cf. (3.39)) if we explicitly choose (cf. (3.17b) for the dehnition of J„) 


£ (n, m, a') 
£ (n, m, a') 


8Nq 4“ I I — 2 [2 (?r 4“ m 4“ u ) J 
8Nq 4“ ItT^Ql 4“3 — 4(?r4“'nr4“ u) 


for a ^ ’ 

for a G > 


(3.43a) 


(n, m, a) ^ 



k 8 iVo 4- I^TVol ~ 5. 


(3.43b) 

(3.43c) 


This choice is also compatible with the “nonlinear” conditions (4.1), which are necessary for the 
treatment of the nonlinearity N'{u) in Section 4. 


3.4. Properties of the parabolic norms and definition of function spaces. In this subsection 
we summarize some of the properties of the parabolic norms lll-IH, HHIIq, and lll-IH^ (cf. (1.16), (3.38), 
(3.40)). 

Lemma 3.4. For given Nq E N and locally integrable u, f : (0, cxd)^ —)■ M such that the generalized 
power series (2.4) exists to order the following estimates (with constants independent of f 
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and u) hold true: 


t 


2i-\-2m—l 




it) 


sup t 

t>0 


2i+2m 


dt™ 
dTu 


^t < |||w|||^ for i G KMo-m \ {0} and m G No, 




it) 


< 


m 


t 


2i-\-2m—l 






it) 


dt < 


for i G and m G No, 

I for ie K^o-m \ {0, /?} and m G Nq. 

(0) 


(3.44a) 

(3.44b) 

(3.44c) 


Furthermore, for any loeally integrable : (0, cxd) —)■ M such that Uq = limx\o (x) exists, we 


have 


u, 


( 0 ) 


< (where the constant is independent ofu^^^). 


Proof. The estimate for has already been proven in [11, Lem. 4.3 (a)]. Estimates (3.44) for the 
coefficients follow quite elementarily by the same reasoning as in the proof of [11, Lem. 4.3 (b)]: 

For estimate (3.44a) we take (a, N) E A with a + iV — m — l = i and i ^ No and obtain®: 

|2 

dx 


d^Ui 

to 

A 

d^Ui 

dt™ 


dt™ 


< 


< 




T. 




j <a-\-N—m—1—5 


dp 


-x^ 


+ 






d^u, . d™n, i 
- -x^ -^- X 




dP 


dr‘ 


dx 


arn 


E 




d^Uj 

_ ^ /J»J 

dt"^ 


dTu 


E 


d^u 


£(A^—m,m,Q;')+4,Q;H-A^—m—1—(5 
2 




dP 


d-x^ 


m,m,Q:')H-4,Q;+7V—m—1+5 


Multiplying with the time weight t2(Q+v)-3 integrating in time, we obtain (3.44a) (cf. (3.38)). 
For proving estimate (3.44b), we take {a, N) G A with a + iV — m — | = i and the same reasoning 
as above (taking the sup in time instead of integrating) leads to estimate (3.44b). The proof of 
(3.44c) is the same as for (3.44a). □ 


Estimates (3.44) will turn out to be relevant for estimating the nonlinearity. However, they are also 
convenient in order to dehne appropriate spaces for our solution and the right-hand side. As in the 
case of standard Sobolev spaces there are two possible approaches: 

On the one hand, one may dehne for a locally integrable fnnction u : (0, cxd)^ -P R its distribntional 
derivatives and take the inhmum over all possible coefficients m* : (0, cx)) —)■ R in the dehnition of 
the norm |||n||| in (3.38). It is apparent from (3.38) that if |||n||| < oo, the coefficients Ui are nniqnely 
dehned almost everywhere. This corresponds to the dehnition of the hL-scale in the standard theory 
of Sobolev spaces. 

On the other hand, it was shown in [11, Lem. B.4] for Nq = 1 that one can approximate any locally 
integrable u : (0, cxd)^ —)■ R with |||n||| < cxd by a seqnence : (0, cxd)^ —)■ R)j^gpj of more regular 
functions fnlhlling: 

(a) G O°°((0, cxd)^) n Oo([0, cxd)^); 


®The assumption z ^ No is merely for notational simplicity, 
the same way for z G Nq. 


The reader may verify that the reasoning works in 
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(b) for every t G [0, cxd) we have + u'p\t)x^ for x <^i, 1, where u^q\u^p'^ : 

(0, cxd) —)■ M are smooth; 

(c) |||m — III —)■ 0 as z/ —)■ cx). 

Taking the closure of all u : (0, cxd)^ —)■ R with (a), (b), and (c) with respect to lll-IH (with Nq = 1), 
we end up with the analogue of what one commonly refers to as the hf-scale of Sobolev spaces. 
Lemma B.4 of reference [11] then states that for iVo = 1 indeed H = W holds®. 

Unlike in [11], where it turned out to be more convenient to rely on the fU-approach of Sobolev 
spaces, we will employ the iL-approach in what follows: 

Definition 3.5. Suppose Nq eN and S > 0 is chosen sufficiently small. Furthermore suppose that 
conditions (3.42) are fulfilled (cf. (3.43) for a specific choice). 

(a) We define the space Uq of initial dataffi^^ as the closure of allu^^l E C°°((0, cxD))nCg([0, cxd)) 

with (t) = Uq^^ + for t 1 with respect to |||■|||Q (cf. (1.16 )). 

(b) The solution space U is defined as the closure with respect to lll-IH (cf. (3.38)j of all u E 

((0 , cxd)^) n Cq ([0, cxd)^) such that u{t,x) = for x 1 , where the Ui : 

(0, CX)) —)■ R are smooth functions of time. 

(c) The space F of right-hand sides f is defined as the closure with respect to lll-IH^ (cf. (3.40) j 
of all f E C°° ((0, cx)2) n Clj ([0, cx)2) such that f{t,x) = EzeA^^uo,/?} fi{t)x^ for X 1, 
where the fi : (0, cx) —)■ R are smooth. 

We mark that in view of Lemma 3.4 the coefficients Wj, Uq^\ and fi of the generalized power series 
are also dehned (at least almost everywhere in time t) for functions u, and /, respectively, for 
which |||m|||, Ill'W^lllg, and |||/|||;^, respectively, is hnite. 

Next we also provide another estimate that in particular guarantees control of the norm supj>o ||m||, 
where 

||m|| := sup |M(a;)| (3.45) 

x >0 

denotes the x-norm in space x. By approximation this also implies continuity of u and derivatives 
if ueU. 

Lemma 3.6. Suppose that Nq E N, m E {0, • • • , iVo — 1}, and i E {O, • • • , £ (l, m, | ± 5) + l} ■ 
Then we have 

supt^™ ||(9 ™D^m||^ + supt^"^ ||(9 ™D^(m — Mo)||^ < |||m|||^ for all u E U, (3.46) 

t>o t>o 

where the constant in the estimate is independent ofu. 

Proof. A proof for an analogous estimate is contained in [11, est. (8.5)] so that we only sketch the 
arguments here. First we may show that, passing to the logarithmic variable s := lux and using 
a cut-off argument in combination with the standard embedding iL^(R) ^ U®(R), the following 
estimates hold 

11^ “ ^o||(o,i] ^ and ||t||[i,oo) ~ Mi-s for any locally integrable v, (3.47) 


®This is also true for the norm IH-Hlo as shown in [11, Lem. B.3]. 
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where the constants are independent of S and ||n||^ := sup^^g^ |n(a;)| for any set A C (0,cxd). 
estimate (3.47) in combination with Lemma 3.4 shows 


supf2™ \\d^D\u-uo) 
t>o 


< 


snp 

t>o 


drD^u-Se,o 


d^uo 

dt^ 


+ supf2"^l|a™D^M 

t>0 


2 

[Loo) 


2 


+ snp 
(0,1] *>0 


d”^no 

dt™ 


(3.44b), (3.47) 

< 


snp 

t>o 



2 

i+l,s 




2 


(3.38) 

< 



2 


By the triangle ineqnality and again using estimate (3.44b), we obtain (3.46). 


□ 


3.5. Rigorous treatment of the linear equation. In this section we prove our main result for 
the linear equation; 

Proposition 3.7. Suppose Nq E N and S > 0 is chosen sujficiently small. Furthermore, suppose 
that conditions (3.42) are fulfilled (cf. (3.43) for explicitly chosen indices). Then for any f E F 
and E Uq there exists exactly one solution u = S f] E U of the linear degenerate-parabolic 
problem (2.1). This solution fulfills the maximal-regularity estimate (3.41). 

of Proposition 3.7. The statement is the generalization of [11, Prop. 7.6] for iVo = 1 to arbitrary 
Nq E N. Since our Banach spaces U and F are nested for increasing Nq, uniqueness is already clear 
and it remains to prove existence. As Proposition 3.7 for Aq = 1 is already proven, there exists a 
unique solution u of (2.1) lying in the space U for Nq = 1. By approximation (cf. Dehnition 3.5), we 
may without loss of generality assume that / G ((0, cxd)^) flCg ([0, cx))) and G ((0, cxd)) fl 
Cq ([0, cxd)), with the expansions u^^^(x) = u® for a; 1, and f{t,x) = \{o/ 3} 

for a; 1, where the /* : (0, cxd) —)■ M are smooth functions of time. By standard parabolic 
theory in the bulk and using [11, Prop. 7.6] also u E ((0, cxd)^) fl Cq ([0,cxd)^) with u{t,x) = 
Uo{t) + ugiffix^ + o (a;^) as a; \ 0. 


A gualitative argument for regularity. Here we argue why the solution indeed has additional spatial 
regularity. As we have observed before, this requires higher regularity in time. At the basis of our 
reasoning is the existence and uniqueness result for the linear problem (2.1) given by [11, Prop. 7.6]. 
In particular we have finiteness of 


,2/3-1 




fc(l,0,/3±<5),/3-l±<5 


dt 


(3.18b) 


,2/3-1 


\dtu\ 


f(l,0,/3±5),/3-l±5 


dt 


and thus we may use hniteness of the norms o/3±5) g-i ±5 some time t = r > 0 to solve 

xdfidtu) + p{D){dtu) = dtf for t > r and a; > 0 

and infer that dtu has additional regularity for times t > r. In principle such a reasoning is possible 
but as it was noted already in [11, Sec. 2], the arguments there^° were also suitable to obtain 


lOyye refer to the discussion in [11, Sec. 2, Eqs. (2.8)-(2.12)] and the (rigorous) time discretization in [11, Prop. 7.6]. 

Utilizing the linear equations (3.5) and (3.6) (where fL) yji^) ami iPe difference between and is immate¬ 
rial due to the choice of the weight 1 -1-5 < 2/3), the integrals in (3.48) are finite through control of the spatial integrals 


r<(i 


1 ( 1)1 


lfc(l,0,l-(5)+4,l-(5 


f^^'^lqi,o.i-i),i-J and/“t(|u;P) 


fe(2,0,5)4-4,1-1-5 


.( 2)1 


lfe(2,0,5),l-|-5 


^ dt, respectively. 
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finiteness of the integrals 




(3.48a) 




k(2,0,S),S 


4 0 *^“0 
t OtU- — 

at 


(3.48b) 


(.(2,0,5),5 
dMQ I 


Now we may use finiteness of the norms g and \dtu —o 5) 5 some time > 0 

with Tjy \ 0 as z/ —)■ cxD. Applying [11, Prop. 7.6] implies that also |||i9tM(- + T;/)||| for A"o = 1 (and 
with appropriately chosen derivative numbers) is hnite^^. In fact, since the arguments in [11, Sec. 7] 
were also applicable for general weights in the interval (|, l), we also obtain 

sup(t -r^)2“-^ 


(t - r,)'“ ' dt < OO 


for all admissible weights a. In particular the integral 






is hnite and we may argue as above that 


(t - r^) - 


\e{ 2 ,i,s),s 


is hnite. Hence ^ _^ < oo and d^u — ^ ^ < cxd for some t = fy > where 

\ 0 as z/ —)■ cxD. Apparently the reasoning can be boot strapped and (taking the limit u —)■ cxd) 


we obtain 


aru-d) 


fc(l,m,a')+2,Q'-i ^ 




, df < cx). 


(3.49a) 


/ |(9rtcd)| , dt < cxD 

for all r > 0, all m = 0, • • • ,Nq — 1, and all admissible a. 

In order to obtain additional spatial regularity, we hrst observe 


(3.49b) 


^ k{2,m,a'),a' ~ 


(3.22) poo 


, dt < cxD for all r > 0. 


Then we may use 


9>(3) ^ p{D - 2)arn;d) 5-^(2) + t, x > 0 


and therefore 


m < om (2) I nofjD')(9""+^7;d) 1 \ff 

t k{2,m,a'),a'-\-l ^ k{2,'m,a^),a'-\-l ^ ^ k{2,m,a'),a^ I * 


im+1 (2)1 


I k(2,m,a'),a' 


(3.18),(3.22) 


k{l ,m+l ,a')+A,a' 


^^Notice that the solution of [11, Prop. 7.6] coincides with dtu, because uniqueness under the assumption 
k 0 i-( 5 ) dt < 00 holds true as can be seen in the uniqueness proof for [11, Prop. 7.6]. 
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Since by construction ,9™7-(2) ^ dt < oo (cf. (3.36)) and by (3.49b), we also have 


Our aim is to show 


1 k(2,m,a'),a'-l-l 

^ ^11 T > 0. 


^ ,,, ,, for t > 0. 

^ /c( 2 ,m,Q:')+ 4 ,Q:' + l ~ ^ /c( 2 ,m,Q:'),Q:' + l 


(3.50) 


This follows by applying elliptic regularity of p{D — 2) (Proposition 3.1), but care has to be taken 
as the decay assumptions (3.23) have to be satished: As in the proof of [11, Lem. B.3] we may 
approximate by a sequence of functions G O“((0,cxd)) in the norm \-\k( 2 ma') a'+v 

Solving p{D — as in [11, Lem. B.3], we hnd a smooth solution such that 

= 0{x‘^) as a; \ 0 and = O as x oo for all i > 0. By Propo¬ 


sition 3.1, 


< 




k(2,m,a'),a'+1 


—)■ 0 as z/, z/' —)■ cxD and hence 


-^L(2,m,a0+4.a'+i ^ 0 as z/ ^ OO for some locally integrable a;. Since L(2,„,„q+4,a'+i ~ 
\^^''^\k{ 2 ,m,a'),a'+v particukr |a;|fc(2,^,„q+4,aMi < I«'),«'+! t > 0 by continuity of the 
norms. Hence, for establishing (3.50) it remains to prove 

u = almost everywhere. (3.51) 

Note that 

p{D - 2) (u- = p{D - 2) lim = lim = 0 

almost everywhere and hence 

(jj — G kerp(Zl — 2) = span ja;^, x^~^‘^, x^~^, x^ |. (3.52) 

Because of L+pji i-5 ^ necessarily = o (a;^”*^) as a; \ 0 and 

a; oo for t > 0. Together with |i^|fc(2,m,o')+4,a'+i ^ therefore 

u — = o as a; \ 0 and u — = o as x oo 

almost everywhere in t > 0. In view of (3.52) this implies (3.51). 

Finiteness of supj>.,- 


, „ ,, 1 for r > 0 can be obtained by the standard trace estimate 

k(2,m,a')+2,a'+^ 


sup 1,(2,^ 

t>T ^ 


,,, ,,, + ) dt. 

^ /c(2,m,Q:'),Q:' ^ /c(2,m,a:')+4,Q:'+l / 


A proof is contained e.g. in [11, Lem. B.2]. 

Apparently the argument can be boot strapped and we obtain 


sup 0^ W k{N-m,m,a')+2,a'-N-m-^ ^ 


t>T 

/*oo 


|om+ly^(A-m) 


k(N—m,m,a'),a'+N—m—2 


and 


Id^W^N-m) 


k(N—m,m,a')+A,a'+N—m—l 


< oo. 

(3.53a) 

dt < oo. 

(3.53b) 

dt < oo 

(3.53c) 


for all r > 0, all A^ = 1, • • • , A"o, m = 0, • • • , A^ — 1, and all admissible a. 
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Proof of estimate (3.33). Now we consider equation (3.16), where and ~ defined 

through (3.18) - are smooth. We are aiming at deriving estimate (3.19). Therefore it is more 
convenient to use the logarithmic variable s := Inx, for which equation (3.16) reads 

^SQ^+I^{n) ^ for t > 0 and s G R. ( 3 . 54 ) 

In order to simplify the notation, we write 

a; := P(C) := p(( - n), and ^ := + e^qn{ds)dT^\^^\ (3.55) 

that is, equation (3.54) can be rephrased as 

e^dtOJ + P{ds)u = (p for t > 0 and s G R. (3.56) 


We take a cut oS rj e C°°(R) with rj > 0, 77 ( 5 ) = 1 for |s| < 1, and 17 ( 5 ) = 0 for |s| > 2 . Let 
Vni-^) •= test equation (3.56) with in L^(Rs), where p is in the coercivity 

range of P(C) (We specifically choose p := a' + n — 1 later on.). Thus we arrive at 


1 d 
2 df 



e-^>^^r]luP{d,)uds 


e-^^^^vlipcods. 


Now we commute one factor rjR with the differential operator P{ds). Since in the commutator of 
rjji and P{ds) at least one derivative acts on rjji, we obtain through integration by parts 


^ '^'^''{riRUjYds + [ e ‘^f^''{T]Ru)P{ds){r]Ru) ds 
2 dt jk Jr 

f e-^^^iVRpf ds+"- [ ds + ^ / + {dsuf + ds, 

2£ Jr 2 Jr R J_2r 


where C* > 0 is independent of R and e > 0 is arbitrary^^. As a next step, we use coercivity of 
P{ds) and obtain 


A / e ^^)"{r]RUjfds + [ e {{r]RUjf + {dsPRuf + (d'^^pRuf) ds 

Jr Jr 


< 

r\j 


e-'^^(w)'ds + A 



[uj^ + idsUiy + idluiy) ds. 


(3.57) 


Estimate (3.57) is a basic estimate that has to be combined with a higher-order estimate to arrive 
at (3.19). Therefore, we again consider equation (3.56) and apply the operator e~‘^^^{ds — IY^'^pr 
to it (where k = k{n,m,a')). Testing with in /.^(R^), we obtain after integrating by 

parts and using a standard interpolation estimate 

A / e<^-"^>{d:^YvRUj)Yds+ [ e-^^^ (df^YvRu))" ds 

Jr Jr 

K p k +4 p -| p2R 

<£-'V / ds + sV / e-2<“(S)(,„w))"ds+- / e""'-* Q" + (SJ+YQ ds, 

Jr j^Q Jr ^ J-2r 

(3.58) 

where e > 0 is arbitrary. By interpolation and increasing C, the second term in the second line 
can be absorbed. Using (3.57), (3.58), and interpolating once more, we obtain, after undoing the 
transformation s = Inx, 

^ \vrYI+2,^-i^ + \vrYI+4,^ Y IwIA + \v2rYI+a,^,- 


12 


All estimates and constants in this part only depend on k and r . 
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Multiplying with where a :=a + n + m — ^ and integrating in time (assuming r' > r > 0), we 
obtain 

sup \VRUj\l+2.u-h + f \VR^\l+4.u^t 


tS[r,T'] 


< 


|hRW|i=rL+2,M-i + / lwL,^dt + 2a 1 dt 


1 

'r 




(3.59) 


Due to (3.53) and (3.55) all terms appearing in (3.59) remain hnite in the limit R ^ oo (note that 
lh 2 Ri^L_L 4 „ 5) IdRw^'^A,, ,,,, ,, , for t > 0) so that we obtain 

I \K+4,fi rsj I t \k{n,m,a')4-4:,a'4-n—l ' 


sup \uj\l^ 1 + / 

te[r,T'] ^ Jt 


a<T \, ,|2 


< 


_2cr 


^\t=T 


k+ 2 ,^—-i ^ 


2a I ,2 

I K,/i 


dt + 2(7 / ^ lea 


2a-l I, .|2 1, 

K+2,/.-i db 


Undoing the notational change (3.55) and increasing the set [t, t'], we end up with (3.19), provided 


that for a + n + m > | 


^ 2(7 


t^lt= 




_ ^2(a+n+?n)—3 


^|i=r 


fc(n,m,a')H-2,Q;'+n—^ 


^ 0 as r \ 0 (3.60) 


holds true at least for a subsequence. This is already clear for (n, m) = (1,0) and a' G (i, l) from 
the extension of the discretization argument in [11, Sec. 7]. From (3.19) for {n,m) = (1,0), we get 


Utilizing 

foo (3.28) r°° 2 / 1 \ 

k^' 1 fe( 2 ,o,«')+ 2 ,a'+i ~ N^'l(iAc.'+|)+ 4 ,a'+| for a'G (^0, - j , 

we conclude that (3.60) holds for {n,m) = (2,0) and a' G (O, |), too. Additionally, we have 


(3.29) 




fc(l,0,«'+i),(a'+i)-l 


for a' G (O, |), that is, (3.60) also holds for {n,m) = (1,1) and a' G (O, |). 
Next, (3.19) yields 


(3.61) 

dt < CXD 


dt<oo for a' G (^0, 0 . (3.62) 


As 


(3.27) 


^2a+l ^ (1) 2 , 1 dt < / 


fc(l,l,«' —|)+4,o'—i 

for a' G (|, l), (3.60) is true for {n,m) = (1,1) and a' G (|, l) and (3.19) yields 

/ ^^<00 for u' G 1 


(3.62) 

dt < OO 


(3.63) 
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Since we also have 


> „ (3.22) roo /o go) 

^2a+l|^ (2)1 , , dt < / , dt < OO 

I ^ I/c(2,0,q:') + |/2| —1,0:^ / ' \k{l,l,a')-\-4,a' 


for a' G (O, |), we get from (3.19) 


1 ^( 2 , 0 ,«'),«' + ^ «' ^ (O’ ^ ' (3-64) 


This enables us to estimate 




fc(2,0,a')+2,a'+l ~ 


(3.27) /‘°o . , X 

, , dt < / t2(«-2)+l |y,(2)| 


fc(2,0,Q'-5)+4,(Q'-L)+l 


(3.64) 

dt < CX) 


for a' G (4, l), thus proving (3.60) with {n,m) = (2,0) and a' G (4, l). 

Induction on n + m verihes (3.60) for all a + n + m > |. 

Now, starting from (3.19), we follow the reasoning of Section 3.3 and arrive at estimate (3.33). 


Proof of estimate (3.41). To obtain the desired maximal-regularity estimate (3.41), we only need 
to be able to apply Proposition 3.1, that is, we need to verify (3.23). We exemplify the reasoning 
by treating the last line of (3.35). For a hxed value (a, A^) G ^ it suffices to show 

D^df^u = ^ i^P' + o {x^) as t \ 0, (3.65) 

i<e 

where := a' -|- iV — m — 1, a' = a ± 5 G (0,1), and £ = 0, • • • , /c (iV — m, m, a') -|- 4. 

Due to |||m|||^ < OO, we have = o{x^) as x \ 0. This implies D^d'f-u{x) = 

D^ 7 r(”^)(a;) -|- D^u{x), where D^u{x) = o(x^) as x \ 0 and 7 r^”^)(a;) is a solution of the homoge¬ 
neous equation P{D)n^^'> = 0 with P(C) ;= (nf=o”'"V(C -/i^)) (nieJ^_„(C - *)) (cf- (3.18b)). 

Hence , where the Uj are smooth functions of time t and j is a zero of the polyno¬ 

mial P{C)- In fact, the set of admissible exponents j is smaller: Since we know that Proposition 3.7 
also holds for iVo = 1, we need to have 7 r^”^^(a;) = + ^^^x^ -f o{x^) as t \ 0. Hence j G {0, /3} 

oi (3 < j < Nq. 

Now suppose that /3 < j < Nq and j ^ K^o- Since P does not appear in the expansion of /, 
we obtain from the linear equation ( 2 . 1 a) + Pij)uj = 0 , that is, if ^ 0 necessarily 

^ ^ 0. This implies j — [jj G {0,/9} (as otherwise negative j would appear) thus contradicting 

the assumption j ^ iPjVo (cf- (2-3)). 

Hence indeed (3.65) holds true and Proposition 3.1 yields estimate (3.41). □ 


4. The nonlinear problem 

4.1. The main results. In this section we prove our main theorem. For the treatment of the 
nonlinear problem (1-11), we need additional assumptions on the numbers i{n,m,a) which we 



30 


MANUEL V. GNANN 


already state here: 

i{N — m, m', a') > i {N — m,m, a') for m' < m — 1, 




i {N — m, m, a') 


+1 where 


0 < m' < m, 

a' + N — m — l</3, 

0 < m' < m, 

a' + N — m — l>/3, 


i ^1, m', ^ > i {N — m, m, a') + 3 where 

i — m, m, a' — > ^{N — m, m, a') + 2 for all a' G , 

i — m — 1, m, a' + -^ > i {N — m,m, a') + 2 for all a' G ^ 0 , , 


k > i{l, m, a') — 2 for all a' > (3 and m > 0 , 
£( 1 , 0 ,1-5) > £(l,m,a')+ 2 ,' 

£(2,0,5) > £(l,m,a') + 2 


for all a' > (3 and m>l. 


(4.1a) 

(4.1b) 

(4.1c) 

(4.1d) 

(4.1e) 

(4.1f) 

(4.1g) 


For the explicit choice (3.43), we are indeed able to fulhll the “linear” conditions (3.42) and the 
“nonlinear” conditions (4.1). The relevance of conditions (4.1) will become clear in the proof of 
the main estimate for the nonlinearity (cf. Proposition 4.2). The aim of this section is to show the 
following statement: 


Theorem 4.1. Suppose Nq G Nq and further suppose that conditions (3.42) and (4.1) are fulfilled 
(cf. (3.43) for an explicit choice). Then there exist e > 0 and 5 > 0 such that for any G Uq 
with 111'^*'°^IIIq < £ (cf. (1.16)j problem (1.11) has a unique solution u E U. This solution obeys the 
a-priori estimate |||m||| < ||l'*^®|||o ('^f- (3.38)j with a constant independent of 

It is quite apparent that Theorem 4.1 (with Nq replaced by iVo-|-l) yields Theorem 2.1 (the regularity 
and decay properties of = ufit) and Rjsiq = Rjsifixfi) in (2.4) follow from Lemma 3.4 and the 
dehnition of lll-IH, cf. (3.38), by approximation with smooth functions as given in Dehnition 3.5). 
Theorem 4.1 in turn has already been proven for iVo = 1 in [11, Th. 3.1] so that in particular 
uniqueness holds true. The fact that existence also holds in a smaller space follows by applying the 
following proposition: 


Proposition 4.2. Suppose that Nq eN and 6 > 0 is chosen sufficiently small. Furthermore suppose 
that conditions (3.42) and (4.1) are fulfilled (cf. (3.43) for an explicit choice). Then the following 
estimate for the nonlinearity Af(u) (defined in (1.13) and (1.9)j holds true: 

lll■^('*^)llll ~ u ^ U, (4.2) 

m=2,5 

where the constant in the estimate is independent ofu. 


Before proving Proposition 4.2, we show how it can be used to prove Theorem 4.1: 

of Theorem 4-1- In Proposition 3.7 we have constructed a solution operator S : Uq x F ^ U to the 
linear problem ( 2 . 1 ), so that the nonlinear problem ( 1 . 11 ) is equivalent to 

M = 5 [m®,A£(m)] . (4.3) 

This hxed-point problem can be solved by applying the contraction-mapping theorem, which for 
TVo = 1 has indeed been carried out in the proof of [11, Th. 3.1] under the assumption of small data 
This implies that for |||'U*'°^|||p 1 the sequence dehned through 

:= 5'[m®, 0] and := S [m^°\ A£ for z/ > 1 


(4.4) 
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converges in U for Nq = 1 to the unique solution u. In particular converges point-wise to u 
(cf. Lemma 3.6). For arbitrary iVo > 1 we note that since G U, applying the maximal-regularity 
estimate (3.41) (cf. Proposition 3.7) and the nonlinear estimate (4.2) (cf. Proposition 4.2), we obtain 


llh^'^lll <C'|||m(°)|||q and |||m("+i) 111 <C'(^|||M(°)|||p + ^axJl|^ for z/> 1, (4.5) 

where C > 0 is independent of Assuming e < (4C^)“^ and 2C > 1, estimates (4.5) upgrade to 
||h^'^|||<^ and < ^ + C'^axJ||M(")|||" for z/> 1. (4.6) 


Inequalities (4.6) in turn imply 


u 


for all z/ > 1. 

ZLy 


(4.7) 


By weak-=t:-compactness, has a subsequence that converges in the weak-=t:-topology of U to the 
unique limit u and therefore u E U. By weak-* lower semi-continuity of the norm lll-IH, estimate (4.7) 
implies 


lull! < 


(4.8) 


Applying the maximal-regularity estimate (3.41) (cf. Proposition 3.7) and the nonlinear estimate 
(4.2) (cf. Proposition 4.2) once more to the hxed-point equation (4.3), we obtain 


u\\\<C 



0 


-|- max 

m=2,5 





which, in view of (4.8) and since 2C > 1, upgrades to |||m||| < 


□ 


4.2. Nonlinear estimates. We conclude the paper with estimating the nonlinear part J\f{u) of 
equation (2.1a): 

of Proposition 4-2. We repeat some of the observations on the formal structure of the nonlinearity 
Af{u) that are contained in [11, Sec. 8]: 


The formal structure of the nonlinearity. We recall that the nonlinearity Af{u) can be written as 

J\f{u) p{D)u - A4sym(l + «,••■ A + u), 

where p{D)u = 5A4sym(l + uAr '' A) and Aisym is the symmetrization of Ai. Since we know^^ 

A4,y^(l, • • • , 1) 0, (4.9) 

we can conclude that by multi-linearity M{u) is a linear combination of terms of the form 

A4sym (m,m, with G {1 ,m}. (4.10) 

As u can be decomposed into u = {u — uo) -l- uq, once more appealing to (4.9) and using multi- 
linearity, we infer from (4.10) that J\f{u) is a linear combination (with constant coefficients) of terms 
of the form 

Aisyra {u — with E {uq , U — Uq} , E {I, Uq , U — Uq} . (4.11) 


^^The identity (4.9) holds, as the single factor D appears once, because the traveling wave utw = 0 is a solution 
of the nonlinear problem (1.11). 
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For the first entry of A^sym in (4-11) we may further use u — Uq = {u — Uq — upx^) + upx^ and the 

fact that A4sym(a^^, 1, • ‘ ‘ > 1) = \p{D)x^ 0. Thus A/'(m) is a linear combination (with constant 

coefficients) of terms of the form 

-Alsym {u — Uq — U^X^ , • • • , with G {Uq, U — Mq}, G {1, Mo, U — Mq}- 

(4.12a) 


or 


■Msym {u/SX^,u — with E {l,Uo,U — Uq}. (4.12b) 

The norms. In view of the dehnition of the norm lll-IH^ in (3.40), it suffices to estimate the expression 


^2(a+A)-3 


d^Afi 


u 


E 

^<i<a^+A^—m —1 


d"(V(«.)). , 
di™ 


dt. 


(4.13) 


£,a'+N—m—l 


where i = i{N — m, m, a') and as in (2.4) 


Af{u) = {Af{u))-x^ + O [x^°) as a; \ 0 and t > 0, 

*SAjVq 

for each (a, TV) G A, with a' = a ± <5 G (0,1) and m E {0, • • • , iV — 1} individually. We distinguish 
between sub- and super-critical terms, that is, terms for which either a'-l-iV — m — 1 < /Tor 
a' + N — m — l>/3. 

Estimating the sub-critieal terms. In this part of the proof, we concentrate on the suh-critieal 
terms, that is, the terms for which a' + N — m — l<j3. Then the sum 

is identically zero. It is convenient to use the decomposition (4.11) for the nonlinearity A/'(m), so 
that by distributing the H-derivatives and applying the triangle inequality, it suffices to estimate 


^2(a+A)-3 Qrn / (1) . . .^(5)^ 2 ^ 

tv / a.' -\-N—m—1 ’ 


where 


u 


= D^^{u-uo) with 0 < < £ + 4, 


G {mq, {u — Uq) } with 0 < £2 < 


£ + 4 


G ^l,Uo,E^^'(u — Mo)} with 0 < £j < 


£ + 4 


for j > 3. 


Now we distribute the time derivatives on the individual factors, so that it suffices to estimate 

r°°^2(a+A)-3l (1)...^(5)|2 , 

/ I la'+A'—m—1 ’ 

^0 

with new factors obeying 

(^(1) ^ dp^£)^^(u — Mo) with 0 < £1 < £ + 4 and 0 < mi < TV — 1, 


LO 


f d™27/„ 1 / -L 4 

^ I A rn 2 ’ '9r^-D^^(M — Mo) } with t) < £2 < —-— and 0 < m 2 < TV — 1, 

I L16 J ^ 


e { 5^.0, 


d”^^Mo 

dt^i 


, (m — Mo) } with 0 < £j < 


£ + 4 


and 0 < mj < iV — 1 for j > 3, 
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where ~ This enables us to obtain the following bound: 


^2(a+A)-3l (l)...^(5) 


^ dt < / t2(a+A-m+mi)-3 ^(1) 

o'+A—m—1 ~ / 


o'+A—m —1 


sup t 




(4.14) 


where ||n|| := sup^-^g |u(t)|. Then we note that 

^oo pco 

/ .2(o+A-m+mi)-3 I (1) |2 r]f< j.2(«+A-m+mi)-3 I ami |2 i, 

^ ^ ^ “h(A-m,mi,a')+4,a'+A-m-l 

Jo Jo 

^ III ll|2 

for the first term on the right-hand side of (4.14), provided that condition (4.1a) is fulhlled. 

We now turn our attention to the terms with j > 2 in (4.14). If we have 

suPt>o^^”^^ ||a;*^T||^ < |||m|||^ by estimate (3.44b) of Lemma 3.4. Finally, if D^i{u — Mq), 

Lemma 3.6 yields supj>gt^™'J' ||^ < |||'*^|||^; provided that condition (4.1b) is fulhlled. 

In summary we can estimate the term in (4.13) in the sub-critical case a'-l-iV — m — 1 < /Sas 
follows: 


^2(«+A)-3 




dt<iii«r X (i+iii«r), (4.15) 


I3<i<a'+N—m—l 




where £ = — m, m, a'). 


Estimating the super-critical terms. We now aim at estimating terms of the form (4.13) with ^ = 
i{N — m, m, a') for each (a, N) E A with a' = a ± 6 E (0,1) and m E {0, • • • , iV — 1}, where we 
assume a'-l-iV —m —1 > (3 (the “super-critical” terms). Here we use the decomposition oiJ\f{u) in 
the form (4.12), that is, the super-critical terms (4.13) can be estimated by a sum of terms of the 
form 


^2(a+A)-3 


d"‘ 

A+m 




df, (4.16) 


i<a'+N—m—1 


' Oi'-\-N—m—1 


where one of the following two situations occurs: 

(a) The obey (cf. (4.12a)) 

^(1) ^ i^u-UQ- ugx^) with 0 < < £ + 4, 

e {mq, {u — Uq) } with 0 < £2 < ^ + 4, 

g |1 ,mo,Z1^^(m — Mq)} with 0 < f'j < — for j > 3. 

(b) The caOl obey (cf. (4.12b)) 

a;d) = ugx^ 


;(2) = £)^ 2 (m_Mo) with 0 < £2 < ^ + 4, 


£ + 4 


G { 1 , Mq, (m — Mq)} with 0 < £^ < —^— for j > 3. 
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Estimating (4.16) for terms of the form (a). As a first step, we distribute the time derivatives on 
the individual factors so that terms of the form (4.16) fulhlling (a) can be estimated by a sum 
of terms of the form 


( 2 („+K)-S ^( 1 ) . , . ^(5) _ 


(4.17) 


i<a' -\-N—m—l 




where the are given by 

(^(1) ^ — Uq — ugx^) with 0 < < £ + 4 and 0 < mi < m, 

r 1 

G < — -, — uo) > with 0 < £2 < f' + 4 and 0 < m 2 < m, 


d^^uo 


G h 


d”^%o 


, D^^{u — no)l with 0 < ij < ^ ^ and 0 < m^ < m for j > 3, 


with Since /9 > | we have = 0 for i < 1 (cf. (2.3) for the dehnition of admissible 

exponents i). Hence we may decompose according to 


cn«= cn«- 




l<i\<a.' -\-N—m—l 






l<ii<a.'+N—m—l 


and instead of terms of the form (4.17), we may estimate terms of the structure 


^2(«+V)-3 I ^( 1 ) 


l<i\< ol' -\-N—m—l 




' a'-\-N—m—l 


(4.18a) 


‘^2(a+A)-3 ^(1) ' ^(2) . . .^(5) _ ^ (^(2) . . (4_18b) 


i<a.'+N—m—l—i\ 


' 0 !.'-\-N—m—l—ii 


respectively, where in the latter case 1 < ii < a' + N — m — 1. Furthermore, the term (4.18b) 
identically vanishes if G [hmjfl, all j > 2, so that we may assume — 


Uq) in this case. 

We begin by estimating terms of the form (4.18a) through 


^2(a+JV)-3 I ^(1) 


l<ii<a.' -\-N—m—l 


a)^il \ 7;(2)..../5) 


X " I o;' ' • • * o;' 


'a'-\-N—m—l 


^ / ^2(a+7V-m+mi)-3 ^(1) 




dt X JJ^supt 


lh,(i)| 


l<ii<a.'+N—m—l 


'a'-\-N—m—l 


j=2 


Then we note that 


^2{a+N-m+mi)-3 ^( 1 ) _ 




l<ii<a'-\-N—m—l 


^2(q:+A^— m+mi)—3 




'a'-\-N—m—l 

|2 


dt < |||m||| , 


ii<a' -\-N—m—l 


' i-\-A,a'-\-N—m—l 
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provided that the indices meet condition (4.1a). By Lemma 3.4 and Lemma 3.6, respectively, we 
can further bound j > 2 if ^ 1 and we fulhll condition (4.1c). 

Altogether, terms of the form (4.18a) can be bound as follows: 


^2{a+N)-3 


U 


( 1 ) _ 


E 




l<ii <a'-\-N—m—l 


dt < lll^lf X (1+ \\\u\f) 


a'+N—m—l 


(4.19) 


For estimating the term (4.18b), we use the LAbound on i.e.. 


roo 

/ ^2(o+A^)-3 

d;> 

2 

LA... LA- ^ 

/o 



i<a' -\-N—m—l—ii 


dt 




poo 

1 j.2mi+22i —1 


Jo 

n 


dt 


(4.20) 


X supfA“+^-"^i-*i)-2 
t >0 


a;(2)...a;(5)- 






i<a' -\-N—m—l—ii 


a' -\-N—m—l—ii 


Then we observe that by estimate (3.44a) of Lemma 3.4 we have for the second line in (4.20) 


t 


2 mi+ 22 i — 1 


UJ. 


( 1 ) 

h 


dt= t 


2 mi+22i — 1 






dt < IIImII 


due to A + mi < iV < Nq. Furthermore, we may use the decomposition 


^( 2 ) ^ ^( 2 ) _ .,( 2)^*2 


E 

^<i 2 <<y' -\-N—m—l—ii 


W-'I" 1 + 


E 

/ 3 <i 2 <a'+N—m—l—ii 


( 2 ) 12 
^2 


and therefore the third line in (4.20) reduces to estimating 


sup 

t>o 


and 


^2(a+A^—^i)—2 


,( 2 )_ 


o;- - 




0 <i 2 <(^' -\-N—m—l—ii 


X 


a' -\-N—m—l—ii 


o 

nsupILL)II" (4.21a) 
t>o 


i=3 


supt 

t>0 


2m2+2i2 


2 

U- 

*2 


X sup 

i>0 


^ 2 {a+N—mi —m2—11—12)—' 2 , 




E (' 




). 


i<a.' -\-N —22 


* a.'-\-N —22 

(4.21b) 

where the term (4.21b) identically vanishes unless caD = D^^{u — Uq) for at least one j > 3. We 

will assume without loss of generality o;)^) = — uq) in this case. Then we can decompose 

a;)^) into 


la 


E 

l3<is<a'-\-N—m—l—ii —22 


df I + 


E 

/3<*3<a'+V—m—1—ii—*2 


(3) 13 


and argue as in the previous step. Hence, in order to close the argument, we need to estimate terms 
of the from (4.21a). Note that not necessarily a + N — m — \—ii G iFiVo but we do know that there 
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exist ivj , ^ 0 : 2 , n'^ G -4, such that 


ai + N — 1712 — - <a + N — m — ii<a 2 + N — m 2 — - 


(4.22) 


and aj + iV is maximal {j = 1) or minimal {j = 2) with this property. Since ii > 1, in particular 
aj < a + ^ and < iV — 1 or a^ < a — | and N < N for j = 1, 2. Further noting that 




we can estimate 


^^2{a+N—m+m2—h)—2 ^{2) _ 


( 2 ) i2 

ul’x^ 


0<i2<ct' -\-N—m—l—ii 


'a'+N —m—l—ii 


^^2(a+N-m+m2-ii)-2 




i2<(y'+N—m—l—ii 


'i-\- 4 ,a'+N —m—l—ii 


(3.18b),(3.39) ^ 

< sup t2ia+N-m+m2-h)-2 I 

^ t>0 


k{N—m,m,a')+ 4 ,a'+N—m—l—ii 


(4.1d),(4.1e) , - ,2 

t>0 kyN—m2,'m2,Ci[j+2,a[+N—m2 — ^ 


+ supt 2 ("^+^)"^ 07^^ 


\k(N—m2,m2,aL)-\-2,aL4-N—m2 — 7 


(3.18b),(3.39) , 

< supt"(“i+'^)-" dr^u 

t>0 




*2<«'i+N— m2 —I 


h ( A—m2 ,m2 ,o^ ) +2,a'j + A—m2 — I 


supd(“+U-= W ^ 

t>o ^ 

i2<Oi2+N—m2 — ^ 


' £(^N—m2,m2,a2)-\-2,a2~^N—m2 — ^ 


where Proposition 3.1 (elliptic maximal regularity) as well as conditions (4.1a), (4.Id), and (4.1e) 
have been used. 

Finally, the product 11^=3 11'^*''^^iT (4-21) can be estimated as before by using esti¬ 

mate (3.44b) of Lemma 3.4 and Lemma 3.6, respectively, provided condition (4.1c) is satished. 
This shows that we can bound terms of the form (4.18b) as 


‘ ^2(a+A)-3 ^(1) ' ^(2) . . . ^(5) _ ^ (^(2) . . . ^(5)^^ 


i<a'+N—m—l—ii 


' a'-\-N—m—l—ii 


(4.23) 


< lllulf X (l-f- lllulf) . 

The combination of (4.19) and (4.23) gives 


(2(„+IV)-S 


di<IIMIl‘x(l + |||u|f) (4.24) 
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for super-critical terms of type (4.17). 

Estimating (4.16) for terms of the form (b). Distributing the time derivatives on the factors 
in (4.16) meeting (b), it suffices to estimate terms of the form 

|2 

(4.25) 

ol'+N— m—l—l3 

where the obey 

A^l 7 / 

= —r—— with 0 < rui < m, 


poo 

/ t2(a+Y)-3|^(l)|2 


Jo 

i<a'+N—m—l—l3 




u 


( 2 ) ^ — Uq) with 0 < £2 < i + 4: and 0 < m 2 < m, 


e <! h 


^rrij ,0? 


d”^^uo 

dt^^ 


, D^^{u — uo) (■ with 0 < £j < 


£ + 4 


and 0 < mj < m for j > 3, 


where Y^j=i '^j ~ = 'nr = iV — l = iVo — 1, we may use the L^-bound on and obtain 

for the term (4.25): 


poo 

/ t2(a+Y)-3|^(l)|2 

1^(2)...^(5) _ ^ 

Jo 

i<a' -\-N—7n—l—^ 


dt 


< 

rs_/ 


^2(/3+A)-3 |c^(i)Mdt xsupt^D-^ |a;(2) 

t>0 


a'—f) 


X 


i=3 """ 


\u 


U) 


a'+N—m—l—g 

ii 2 


(4.26) 


t>0 


Then we note that we can estimate the first term in the second line of (4.26) through 


^2(/3+V)-3 1 (l)|2^^^ / t^(h+N)-3 


dY -1 


Ug 


dt ^-1 


dt < lllul 


by estimate (3.44a) of Lemma 3.4. Furthermore, since we are in the super-critical case, we have 
{t/xY^°‘~^'> <1 + t/x and we can bound as follows^"^ 


sup< supi"l“-'’> |B(t, - t,„)| 


t>0 


i>0 


< 

r\j 


£+3,0:^—/3 
2 


sup |m |^+4 + sup \u - Uolf+ 4,5 + (1 - <^l,v) supt |m - Mo|f+ 4 ,i ±5 


< 


u 


t>0 


t>0 


t>0 


where Proposition 3.1 has been used and conditions (4.If) and (4.1g) need to be satished. 


The product 11^=3 s^Pt>o (4.26) can be estimated by employing estimate (3.44b) of 

Lemma 3.4 and Lemma 3.6, respectively, if condition (4.1c) holds true. In summary we obtain 


r°°^2(«+A)-3|^(l)|2 


Jo 



a'+N—m—l—/3 


(4.27) 


< 


iuif X (i+iiiuiir) 


in the case mi = m = N — l = No — 1. 


the special case = 1 (i.e., Nq = 1 by assumption) we have a = P and the second summand on the right-hand 
side can be dropped. 
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In all other cases, we can estimate (4.25) by taking the C'^-bound on and get: 

|2 


a'+N—m—1—(3 


/'“t2(a+A)-3|^(l)|2 


/O 

i<OL'+N—m—l—j3 




^ sup 

t>0 


(4.28) 


POO 

j ^2{a+N—mi—/3)—3 


/o 

i<a'+N—m—l—l3 


dt. 


a'+A—m—1—/9 


Then we notice 


sup 

t>o 


^2m,+2P ^(1) 2 = 


t>0 




dt"^i 


< 




by estimate (3.44b) of Lemma 3.4. For the last line of (4.28) we may use the decomposition 




so that it suffices to estimate 




I + 


l3<-i2<(y.' -\-N—m—1—p 




j2{ct+N—m+m2 —/?)—3 


UJ 


( 2 ) 


E 

I3<i2<<y.'+N—m—l—(3 


^2 


dfxj^supt^™^ ||cn(^')|f (4.29a) 

o'+A—m—1—/3 i—3 


t>0 


and 


poo 

1 ^27712+2^2 — 1 

(2) 

/o 

12 


X sup 

i>0 


dt 


j2{a+N—m\ —m 2 —/3—* 2 )—2 




E 

i<a.' -\-N—m—l—l3—i2 




* a'-\-N—m—l—(3—i2 

(4.29b) 

where the term (4.29b) is identically zero unless D^^{u — Uq) for at least one j > 3. Since 

3 + ^2 ^ 1, the expression (4.29b) can be treated in the same way as the second and third line of 
(4.20). 

For the term (4.29a) we can treat the product 11^=3 supt>o ||a;^t) ||^ ag before (by using Lemma 3.6 
and estimate (3.44b) of Lemma 3.4, respectively) given that we fulhll condition (4.1c). For treating 
the L^-part in (4.29a) we pick ("Q;i,iv') , (a 2 ,N') E A such that 


Pi + iV — 7712 <a + N — m — 3 < a 2 + N — m 2 (4.30) 

and aj + iV is maximal {j = 1) or minimal {j = 2 ) with this property. Since /9 > | we have in 
particular aj < a + ^ and N — m 2 < N — m — 1 or aj < a — ^ and N — m 2 < N — m for j = 1,2 
under the assumption of a sufficiently small 5 > 0. Utilizing 
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we obtain the estimate 


poo 

1 ^2(q: + TV—777+7712—/3)—3 

a;(2) _ ^ 

Jo 

/3 < 72 < a ^ + TV — 777 — 1 —/3 


dt 


a'+N—m—l—j3 
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^2(a+N—m+m2—0)—3 




E 


(3.18b),(3.39) roo 

< / 


^2(a+A-m+m2-/3)-3 \Qrn2^(N-m) 


i 2 <ct'+N—m—l—P 
2 




dt 


(4.31) 
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^2(ai+A)-3 




N—m2 ] 


fc(A—m,m,a')+4,o'+A—m—1—/3 

dt 


(.+4:,a'+N—m—l—(3 

dt 


/ ^2(a2+Y)-3 

Jo 

(3.18b),(3.39) roo 


o^ 2 ^(A-m 2 


kyN—m2 ,m2 ,o^ j +4,Oj + A—m 2 — 1 
2 
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^2(ai+A)-3 
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/c ^ 7 ^— 7712 ,^ 2 , 0 : 2 )+4,0^2+— l 
2 


dt 


*2<Q;i+A— m2 —1 


dt"^2 


— 


dt 


/ ^2(Q:2+iV)-3 


^ ^ TV — 7712 ,7712,0 ]^ )+4,+ TV — 7712 — 1 
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ar% - 5^ 


i2<a2+N—m2 — l 


dt "^2 


dt 


£( A—m2,m2,«2)+^i‘^2A^~™2 —1 


(3.38) 

< 


U 


where we have used elliptic maximal regularity in the form of Proposition 3.1, (4.1a), and the 
conditions 

N — m,m,a -] > i {N — m,m,a') for all a E ( -, 1 ] , 


(4.31a) 

£ — m — l,m,a + ^ > £ {N — m,m, a') for all a E ^0, . (4.31b) 

These conditions in particular hold if the stronger conditions (4. Id) and (4.1e) are fulhlled. 

In summary we obtain (4.27) for the case of arbitrary mi and N as well. 

Conclusion. Gathering inequalities (4.15), (4.24), and (4.27), we obtain inequality (4.2), thus £n- 
ishing the proof of Proposition 4.2. □ 
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